The coherent scattering function in the reptation modeh 
analysis beyond asymptotic hmits 
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We calculate the coherent dynamical scattering function Sc{q, t; N) of a flexible chain of length A^, 
diffusing through an ordered background of topological obstacles. As an instructive generalization, 
we also calculate the scattering function Sc{q, t; M, N) for the central piece of length M < N oi the 
chain. Using the full reptation model, we treat global creep, tube length fluctuations, and internal 
relaxation within a consistent and unified approach. Our theory concentrates on the universal 
aspects of reptational motion, and our results in all details show excellent agreement with our 
simulations of the Evans-Edwards model, provided we allow for a phenomenological prefactor which 
accounts for non-universal effects of the micro-structure of the Monte Carlo chain, present for short 
times. Previous approaches to the coherent structure function can be analyzed as special limits of 
our theory. First, the effects of internal relaxation can be isolated by studying the limit A'^ — » oo, 
M fixed. The results do not support the model of a 'Rouse chain in a tube'. We trace this back 
to the non-equilibrium initial conditions of the latter model. Second, in the limit of long chains 
(M = A/" — > oo) and times large compared to the internal relaxation time {t/N^ — > oo), our theory 
reproduces the results of the primitive chain model. This limiting form applies only to extremely 
long chains, and for chain lengths accessible in practice, effects of, e.g., tube length fluctuations are 
not negligible. 



I. INTRODUCTION 



The equilibrium dynamics of a dense polymer system, 
i.e., a melt, a solution of high concentration, or a free 
chain moving through a gel, is an important topic of poly- 
mer physics. It has been investigated for many years, but 
still is not fully understood. The problem is quite com- 
plex, even if we concentrate on the motion of a single 
chain. Clearly, its motion is strongly hindered by the 
surrounding chains, which the chain considered cannot 
cross. This has led to the idea that the motion of the 
chain is confined to a tube roughly defined by its instan- 
taneous configuration. Thus the tube is assumed to have 
a random walk configuration which changes only by the 
motion of the chain ends. The ends can retract into the 
tube which thus is effectively shortened, and they can 
creep out of the original tube, thus creating a new tube 
segment in some random direction. The interior parts of 
the tube are assumed to be fixed in space until they are 
reached by the diffusive motion of the chain ends. This 
concept of a tube is one basic ingredient of the 'repta- 
tion' 1^,^ scenario, which certainly is valid provided the 
obstacles confining the chain motion form a rigid, time in- 
dependent network. In a realistic system the surrounding 
chains are mobile, which sheds some doubt on the postu- 
lated existence of a well defined tube. Indeed, there exist 
other approaches more in line with standard many 
body theory, which are not based on the tube concept. 

Most work on the reptation model concentrates on 
asymptotic results expected to hold for long chains in 
special time regions (see Sect. II. A). In comparison to 
experiments or simulations, these results often fail on the 



quantitative level j5|, and partly other theories seem to 
be more satisfactory Thus some work |^ has been 

invested to incorporate additional physical effects like re- 
laxation of the surrounding or specific interaction effects 
into asymptotic reptation theory. However, the evalua- 
tion of the pure reptation model outside asymptotic lim- 
its has found only little attention. In recent work 
we presented such a calculation for the motion of indi- 
vidual segments of the chain. We found that asymptotic 
results, which for the quantities considered take the form 
of specific power laws, apply only to surprisingly long 
chains. Large time intervals are covered by crossover 
regions. Our crossover functions compare very well to 
simulations |^ of the pure reptation model, i.e., to the 
motion of a flexible chain through a fixed regular lattice 
of impenetrable obstacles. Furthermore, also results of 
simulations of melts look quite similar to our analytical 
results for short chains. This is consistent with the obser- 
vation that the tube diameter, i.e., the average distance 
among effective obstacles extracted from the simulations, 
is fairly large. Since in the pure reptation model, the tube 
width is of the order of the effective segment size of the 
reptating chain, this implies that to map melt dynam- 
ics on the reptation model, we have to consider a coarse 
grained chain of effective segment number N/Ng. Here N 
is the chain length (polymerization index) of the physical 
macromolecule, whereas the 'entanglement length' A^e is 
the length of a subchain which shows a coil radius of 
the order of the obstacle spacing. In recent simulations 
JlOt of melts, a value N/N^ ~ 14 was reached, far be- 
low the value N/N^ ^ 50 needed according to our theory 
to clearly identify asymptotic power laws. Other recent 
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simulations reach a value N/N^ « 300, but for this 
chain length they cover only times short compared to the 
characteristic time scales of the reptation model. Still, 
the onset of a first power law regime is seen, again con- 
sistent with our theory. Thus, concerning the motion of 
individual chain segments within a melt, there at present 
seems to be no need to invoke other mechanisms than 
pure reptation. 

The motion of specific segments is easily accessible only 
to computer experiments. Physical experiments often 
measure dynamic scattering functions. Asymptotic re- 
sults of the reptation model for the coherent scattering 
function have been worked out previously but 
our analysis of segment motion suggests that an evalua- 
tion of the scattering function outside asymptotic limits 
is needed. This is the topic of the present work. We 
use the same analytical reptation model as in our previ- 
ous analysis We also measured the coherent struc- 
ture function in Monte Carlo simulations, again using the 
same implementation of the model as previously |^ . This 
allows for a comparison among theory and data, where 
all parameters are fixed by our previous work. Some re- 
sults of the simulations will be presented here, but a de- 
tailed comparison of our simulation results to the present 
and previous theories will be presented in a separate, less 
technical paper |l4| . 

In the next section, we briefly review the basic features 
of the reptation model and recall previous results for the 
coherent structure function. In Sect. Ill, we introduce 
our analytical model and outline the structure of our ap- 
proach. In Sect. IV, we consider those contributions to 
the coherent structure function in which the initial tube 
is not yet completely destroyed by the stochastic mo- 
tion. A rigorous analysis is possible as long as end ef- 
fects can be neglected. These end effects, known as 'tube 
renewal' and 'tube length fluctuations', can be treated 
only in some approximation. We here generalize an ap- 
proach which in our previous work gave good results for 
the segment motion. In Sect. V, we compare our rigorous 
results for the motion within the initial tube to those of 
the model of a 'Rouse chain in a coiled tube' Pro- 
nounced differences are found and their origin is clarified. 
In Sect. VI, we derive an integral equation which takes 
complete tube destruction into account. For long chains 
and times large compared to the internal relaxation time 
of the chain, we recover the results of the 'primitive chain' 
model [p"3P], as shown in Sect. VII. Typical numerical re- 
sults of our theory are discussed in Sect. VIII. It is found 
that tube length fluctuations, which have been neglected 
in previous calculations of the coherent scattering func- 
tion, in fact determine the scattering up to times larger 
than the Rouse time. In Sect. VIII, we also present some 
results of our simulations, which compare favorably with 
our theory. Finally, Sect. IX contains a summary and 
conclusions. The full evaluation of the reptation model 
leads to quite involved expressions, and some part of the 
analysis is summarized in appendices. 



II. REVIEW OF THE REPTATION SCENARIO 
AND OF PREVIOUS RESULTS FOR THE 
COHERENT STRUCTURE FUNCTION 

A. Basic dynamics and time scales 

As mentioned in the introduction, the reptation model 
assumes the existence of a tube defined by the instanta- 
neous configuration of the chain together with the sur- 
rounding obstacles. The chain cannot leave the tube side- 
ways since it would have to fold into a double-stranded 
conformation which costs too much entropy. Those parts 
of the chain, which lie stretched in the tube, essentially 
cannot move. In the interior of the tube only little wig- 
gles of 'spared length' are mobile, as illustrated in Fig. 1 
for the special case of a lattice model. These wiggles 
carry out Brownian motion along the chain. If a wiggle 
reaches a chain end, it may decay and prolong the tube 
by its spared length. Chain ends also may produce new 
wiggles which then diffuse into the interior of the tube. 
This shortens the tube by the spared length of the newly 
created wiggle. In the long run, this random motion of 
the chain ends leads to a complete destruction of the ini- 
tial tube. 

This very simple dynamical model involves several time 
scales. It needs a microscopic time Tq until the segment 
motion feels the existence of the constraints. Tq gener- 
ally is identified with the Rouse time of a chain of length 
equal to the entanglement length Nf,: 

Tq^NI . (2.1) 

To is relevant for the short-time dynamics of melts, where 
the tube diameter typically is found to be quite large 
0,111: Ne-^ 10- 40. For reptation, Tq defines the ele- 
mentary time step, since this theory does not deal with 
the unconstrained motion on scale of the tube diameter. 
A second scale T2 is the time a wiggle needs to diffuse 
over the whole chain. Since in the coarse grained descrip- 
tion, the wiggle has to diffuse a distance of N/Ng steps, 
one finds 

T2 ^ To (^^) ~ . (2.2) 

T2 thus is of the order of the Rouse time of the whole 
chain. Finally, the reptation time T3 is needed to de- 
struct the initial tube completely. Reptation theory |^ 
predicts 

as limiting result for long chains. Asymptotically power 
laws as function of t and N are predicted to hold for the 
segment motion or the motion of the center of mass in 
the time windows To < i < T2 , T2 < t < T3 , and 
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Ta <C i . As a typical result we quote the mean squared 
spatial displacement of some bead of the chain: 

f t^/* for To < i < T2 

(r, (t) - r, (0))2 ) ^ <^ (t/Ar)i/2 for « t « Tg 
[ t/N^ for < i . 

Here rj{t) gives the spatial position of bead j at time 
t. The bar indicates the (dynamic) average over the 
stochastic motion and the pointed brackets denote the 
(static) average over all tube configurations. 



B. Previous results for the coherent structure 
function 

We consider a chain of TV + 1 beads {N segments), 
labelled by j = 0, 1, . . . , iV. The coherent structure func- 
tion is defined as 



Sc{q,t,N) 



N 

E 



Jq(r,(i)-rfc(0)) 



(2.4) 



By definition Sc{q,t,N) refers to a single chain. It can 
be measured by appropriately labelling a few chains in 
the system. Reptation results for Sc{q,t,N) previously 
have been derived by Doi and Edwards ||l3|] and by de 
Gennes 

Doi and Edwards have evaluated a simplified version 
of the reptation model, where the internal motion of the 
chain is neglected. The physical chain is replaced by a 
'primitive chain', which only can slide along the tube 
so that all segments experience the same curvilinear dis- 
placement A^(t). This model therefore reduces the dy- 
namics to diffusive motion of the single stochastic vari- 
able . For the coherent structure function, it yields 
the result (see Ref. chapter 6.3.4) 



Sc{q,t,N) 
Sc{q,0,N) 

Sde{Q,t) 



S 



DE 



Q 



2n2 



E 



sm 
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D{Q) (Q2/4 + Q/2 + a 



(2.5) 
, (2.6) 



where i?^ is the radius of gyration and ^ Tj, ^ N^. 
D{Q) is the Debye function: 



DiQ)^^ie-^-l + Q) . (2.7) 
The ap = ap{Q) are the positive solutions of 



Q 



ap tan ctp = — 



(2.8) 



Neglecting all internal motions, the result can be applied 
only for t ^ T2, i.e., in a time regime where the internal 
degrees of freedom are equilibrated. In the limit of large 
wave numbers Q ~ q^R^ I, the result reduces to 



SDEiQ,r) = (2p- l)"'cxp [-(2p- 1)V] 



p=i 



(2.9) 



This is the scattering from that part of the primitive 
chain which at time t still resides in the initial tube . 

The limit q^R^ S> 1 has also been considered by de 
Gennes. Taking the internal relaxation of the chain into 
account, his result for the normalized coherent scat- 
tering function is a sum of two terms: 

SdG{q,t,N) ^ (l-BdG{q)) S^'\q,t) + Bdoiq) S^'Ht,N) 

(2.10) 

BdG{q)^l-^q'Rl . (2.11) 



The 'creep term' S^'^^t, N) is given by Eq. (2.9) and thus 
describes the large time behavior t ^ T^- It tends to 1 
for i/Tg 0. The 'local term' S'^^\q,t) is taken from 
an approximate evaluation of the internal relaxation of 
an infinitely long one-dimensional Rouse chain, folded 
into the three-dimensional random walk configuration of 
a tube of N/Ne segments. The result reads 



S'(^)(q,t) = e*i erfc^/I 



where 



. _ 3 , 2 n2\2 t t 



(2.12) 



(2.13) 



This introduces an additional (/-dependent time scale 



Td 



3 N {q^RlY 



(2.14) 



which in view of R^g ^ TV, and Td ^ is indepen- 
dent of N . Tq governs the relaxation of segment density 
fiuctuations on scale q~^ . In view of ^ 1, Tq is 

much smaller than Td, and for times t Tq, the creep 
term is constant,_ 5('=)(t, iV) « S^''\Q,N) = 1. On top 
of this plateau, S^^\q,t) yields a peak rapidly decreas- 
ing in time. Note that S^^\q,t) for ti ^ I behaves as 
5W(g,t) = {nti)-^/^. The amplitude of the peak is de- 
termined by BdG, which only depends on and Ng. 

Both these approaches neglect end effects like tube 
length fluctuations, which are governed by the time scale 
T2. The approximations involved greatly simplify the 
analysis but are no essential part of the reptation model. 
In the sequel, we present an analysis of the full model, 
accounting for the internal degrees of freedom and the 
finite chain length. Since all the dynamics is driven by 
the diffusion of the spared length as the only stochas- 
tic process, this yields a unified description of local re- 
laxation, global creep and tube length fluctuations. We 
will flnd that tube length fluctuations, in particular, have 
an important influence for intermediate times and chain 
lengths. Internal relaxation, however, is of much less in- 
fluence than the results referred to above suggest. 
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III. FORMULATION OF THE FULL REPTATION 
MODEL 

A. Microscopic dynamics 

We here recall the essential features of our model. A 
more detailed discussion can be found in Ref. The 
chain is modeled as a random walk of N steps of fixed 
length \rj — rj_i| = £o, j = 1, . . . , N . The motion is due 
to the diffusion of wiggles of spared length ig. These are 
represented by particles hopping along the chain from 
bead to bead, with hopping probability p per time step. 
The particles do not interact, and a given particle sees 
the others just as a part of the chain. If a particle passes 
a bead j, it tracks it along by a distance of the spared 
length £s, which is taken to be the same for all particles. 
The end beads j — 0, N of the chain are coupled to large 
reservoirs, which absorb and emit particles at such a rate 
that the equilibrium density po of particles on the chain 
is maintained on average. Keeping track of the change 
of the particle number in these reservoirs, we control the 
motion of the chain ends: creation or decay of a wig- 
gle at a chain end implies emission or absorption of the 
corresponding particle by the reservoir. 

For the motion of beads in the interior of the tube, the 
essential stochastic variable of the model is the number 
n(j, t) of particles which passed over bead j within time 
interval [0,t\. 



(3.1) 



Some useful properties of A-^ ( j, t) are collected in ||^ , Ap- 
pendix A. We also will need the first moment |n(j,i)|, 
which from Ref. §, Eqs. (3.26), (3.27), is found as 



\<ht)\ 



{pMjM'^ [l-Fi(4poAi(j,t))] (3.7) 



2^ 



T 



1 



(3.8) 



where T{y, z) is the incomplete F-function. 

Except for microscopic times t ^ 2/p, can be ap- 
proximated as 



exp 



-4ptsin — 



(3.9) 



so that the theory involves time only in the combination 

i = pt . (3.10) 

In evaluating the theory, we will use t as time vari- 
able. For ■n?{j,t) > 100, which for N > 100 implies 
t > 10**, 7^i(n; j, i) is well represented by a simple Gaus- 
sian 



{n;j,t)^ (2tt n^ij,t)^ 



-1/2 



exp 



(3.11) 



Here n± (j, t) is the number of particles that came from 
the 'left' (/ < j) or from the 'right' (/ > j), respec- 
tively. Consider, for instance, the motion of segment j 
for a time interval in which it stays in the original tube. 
Its displacement in the tube is given by isn{j,t), and 
since the tube has a random walk configuration, its spa- 
tial displacement is given by 



(r,(i)-r,(0))2 )=4 io Hj,t)\ 



(3.2) 



Since the underlying stochastic process is single parti- 
cle hopping, the distribution function of n{j,t) is easily 
calculated, with the result (Ref. B, Eq. (3.22)) 



Pi(n;j,i) = (5„,„(,-4) = e-"'(j'^) /„ ( n^{j,t) ) , (3.3) 

where In{z) is the modified Bessel function of the first 
kind. The second moment {j, t) is found as [see Ref. 
@,Eqs. (3.24), (3.12), (3.9)] 



n2(j,i) = 2pov4i(j,t) 

A,{j,t)^^ + — y (1 

' N 2N ^ ^ 



o^K = I — 4psin — — . 

^ 2N 



(3.4) 



ai) (3.5) 



sm 



2 (TrK_\ 

\2N J 



(3.6) 



Knowledge of 'Pi{n\ j, t) is sufficient as long as we con- 
sider motion inside the initial tube. End effects introduce 
a more complicated quantity. Within time interval [0, i], 
the tube from the end j = is destructed up to bead j< , 
where j< is defined as 



J< 



c(0,t) 



c(0,i), 



max \—n(Q, s)! . 
se[o,t] 



(3.12) 



(3.13) 



Here 



n(0, s) = mo{s) - mo(0). 



where mQ{s) is the occupation number at time s of the 
reservoir at chain end 0. Thus nmax(0,t) is the maximal 
negative fluctuation of the occupation number of reser- 
voir j = in the time interval [0,t\. In (3.12) we also 
introduced 



(3.14) 



measuring all lengths in units of the segment size £o . Sim- 
ilarly, from the other end tube destruction within time t 
proceeds to bead 



(3.15) 
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with n„iax(-^, being the maximal negative fluctuation 
of the occupation number tojv(s) of the reservoir at chain 
end N. The stochastic processes mo(s) or mAr(s) are not 
Markovian, since a particle emitted by a reservoir can be 
reabsorbed by the same reservoir later. This induces a 
correlation which dies out only if the particle has time to 
reach the other reservoir, i.e., on time scale T2. For such 
a correlated process, the distribution and the moments 
of rimax cannot be calculated rigorously, even though ar- 
bitrary moments of n(0, s), involving any number of time 
variables s, can be evaluated (see Sect. III). As soon 
as tube renewal comes into play, we therefore have to 
resort to some approximation. 

Some important quantity entering our theory is the av- 
erage ri,nax(0, t). It, for instance, yields the motion of the 
end-segment via the relation [Ref. Eq. (2.12)] 

( (ro(i)-ro(0))2 ^ = 2£s tl ^WM) • 

We use the expression (Ref. Eq. (5.1)) 

^WM)=E^%^ , (3.16) 

s— 1 

which is correct for a Markov process. Using in Eq. (3.16) 
the exact moments |?^(0, s)\ (Eq. (3.7)), we in essence ap- 
proximate the correlated process by a sequence of Markov 
processes which for each time step s yield the correct in- 
stantaneous value of |7i(0, s)|. This 'mean hopping rate' 
approximation, which was discussed in more detail in [^, 
gives good results for larger times. For microscopic times, 
it underestimates nmax(0,t) by about a factor of 2, but 
with increasing time it approaches the full result for nmax 
as found in simulations. For t T2, the deviation for the 
motion of the end segment, which is most sensitive to our 
approximation, is of the order of 10 % only (see Fig. 9 of 
Ref. @). 

B. Outline of our calculation of the coherent 
structure function 

The basic quantity to be considered, is the scattering 
from a pair of beads 



^((?,t;j,fc,7V) = ^ e^q('^j W-^fc(0)) ^ , (3.17) 

which must be summed over the beads to find the co- 
herent structure function Sc{q-,t; N). To get some infor- 
mation on the contribution of the different parts of the 
chain, we consider a slight generalization in which we 
sum only over the M + 1 central beads 

N + M 
2 

S,iq,t;M,N)^ J2 S{q,t; j, k, N) . (3.18) 



Clearly, the coherent structure function of the full chain 
is 

S,iq,t;N)^S,iq,t;N,N) . (3.19) 

To calculate S{q, t; j, fc, N), we first consider the contri- 
bution 5*^^^ (q, t; j, k, N) which results from those stochas- 
tic motions for which a part of the initial tube still exists 
at time t. (The upper index (T) stands for 'tube'.) We 
then can set up an integral equation for S{q, t; j, k, TV), in 
which S'--'^\q,t; j, k, N) shows up as inhomogeneity (see 
Sect. VI). Furthermore, for t <^ T3, contributions where 
the tube is destroyed completely, are negligible, and S'^'^'^ 
coincides with S. 

iS*^"^^ {q, t; j, k, N) incorporates the effects of internal re- 
laxation and tube length fiuctuations, and its calculation 
is the most tedious part of our analysis. We here need 
to simultaneously control the motion of segment j and 
of the chain ends. More specifically, we will need the 
distribution function 



^maxj("-™>"j;0 = 0(j< -i>) '5n„,nmax(0,t) ^n.Mj.t) , 

(3.20) 

i.e., the simultaneous distribution of n{j, t) and ninax(0, t) 
under the constraint that a part of the initial tube still 
exists. Again the correlated nature of the stochastic mo- 
tion of the chain ends prevents a rigorous evaluation of 
^max j ' random walk theory to construct an 
approximate functional form. The result depends on n,„, 
rij and t only through the rescaled variables n,nl nm{t), 
1/2 

and rij/ {n'j{t)) , and in the spirit of our mean hopping 

rate approximation, we in these variables replace the nor- 

1/2 

malizing factors nm{t) and (n| (t)) of the random walk 

by their counterparts for the proper correlated process. 
In essence, this again amounts to replacing the correlated 
stochastic motion of the chain ends by a whole sequence 
of uncorrelated random walks, parametrized by an ef- 
fective hopping rate p' . This hopping rate is adjusted 
such that the random walk which replaces the correlated 
process for final time t, at that time yields the correct 
moments n„i{t) = 7imax(0,t) and n'^{t) = {n{j,t)y. (It 
in fact yields the correct Gaussian distribution of the 
single variable n{j,t).) As discussed in Ref. Sect. 
V.B, p' changes from a value pop at microscopic times to 
Pqp/N for i ^ T2. Since pop governs the short time mo- 
tion of a segment whereas pop/N is the mobility of the 
primitive chain, the mean hopping rate approximation 
smoothly interpolates between these more rigorously ac- 
cessible limits. This will be discussed again in Sect. IV. C, 
after we presented the details of our approach. 

Our theory involves three important time (and seg- 
ment index) dependent parameter- functions: 

C = C(t) = J^I rimax(0,t) (3.21) 
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measures the extent of tube destruction and thus ac- 
counts for tube length fluctuations. For nniax(0,f), we 
use the approximation (3.16). It turns out that the time 
dependence of c{t) which very slowly tends to its asymp- 
totic limit c{t) const t^/^ (cf. Eq. (7.2)), is responsi- 
ble for the well known crossover behavior of the reptation 
time: T3 ^ iV^"'*, where z^s slowly approaches its asymp- 
totic value ZeS 3 from above. (A detailed discussion 
of the reptation time will be given in a separate paper.) 

A second function, a{j, t), measures the coupling of the 
motion of an interior segment j to the motion of a chain 
end. Initially, this coupling vanishes, but it increases with 
time due to particles created at a chain end and travel- 
ing over segment j. If this coupling is fully developed, all 
segments approximately have moved the same distance 
in the tube and the primitive chain model results. The 
precise definition of a{j,t) is given in Eq. (4.17). 

Finally, it should be noted that the effective mobility 
of a segment for t <^ T2 depends on its position in the 
chain, an effect already present for free Rouse type mo- 
tion. This is taken into account by the function b{j,t), 
which is defined in Eq. (4.36). 

Having described the main ideas of our approach, we 
now turn to the details. We first construct and analyze 
the tube conserving contribution Sc {q,t; M, N) to the 
structure function. 



IV. TUBE CONSERVING CONTRIBUTION TO 
THE STRUCTURE FUNCTION 

In this and the next section, we consider the contribu- 
tion of those stochastic processes, which do not destroy 
the initial tube completely, i.e., we insist on the inequal- 
ity 

.7>-.?<>0 , (4.1) 

since due to the definitions (3.12), (3.15), the tube has 
been destroyed up to segment j< from chain end or 
j> from chain end N, respectively. We first construct a 
formally exact expression for the corresponding contri- 
bution S'-'^\q,t]j,k,N). Its summation over indices j 
and k as in (3.18) and (3.19) yields the tube conserving 
contribution to the coherent structure function. 

A. Exact expression for S^'^^ [q, t; j, k, N) 

Depending on the relation among j, k, j<, and j>, we 
have to distinguish several cases. We use the notation 

j{t)=j + isn{j,t) , (4.2) 

and we illustrate the analysis with two typical cases 
shown in Fig. 2. Fig. 2 a gives a schematic sketch of 
a situation, in which the inequalities 

j{t)<k, j<<k, 3>>3{t) , 

hold. That means that tube renewal from chain end 
has not passed over the original position of segment k, 



and segment j at time t is not found in the part of the 
new tube created from chain end A^. Furthermore the 
new position of segment j, if measured along the tube, 
is closer to the new position of chain end zero than the 
original position of segment k. The relative ordering of k 
and j> , or of j (t) and j< is unimportant. As is clear from 
Fig. 2 a, the path connecting j{t) and k has k—j{t) steps, 
and since the chain configuration is a random walk, we 
find 

^gzq(r,-(t) - rfc(O))^ ^ ^-f{k - j{t))^ (4 3) 

where f = (4.4) 

As a result, the contribution of such configurations to 

S^'^\q,t;j,k,N) reads 

e{k-j{t)) e(fc-j<) Q{j>-m) e(j>-i<) e-^-^c^-^w), 

where the discrete 6-function is defined as 

«(»)^{; : ::"-;:!,:.. 

Now consider a typical case of other type, shown in 
Fig. 2 b. It is defined by the inequalities 

jit) <k , j<>k , 

and differs from the previous one in that tube renewal 
from chain end has passed over the original position of 
segment k. The thus created part of the new tube nec- 
essarily contains the new position of segment j, and the 
random walk connecting j(f) to k has {j< — k) + {j<^—j{t)) 
steps. We thus find the contribution 

e(fc-i(i))e(j<-fc-i)e(i>-i<) e-^(^i<->'-m) . 

The other Cctses compatible with j> ^ are given 
by the relations [ j{t) < k, j> < jX*) ] , [ j{t) > 
k, j< < j{t), j> >k] , [ j{t) > k, j< > j{t) ] , and 
[ j{t) > k, j> < k] . Proceeding as above, we after some 
manipulations with the G-functions arrive at the result 

{q, t;j,k,N)=Si (g, t;j,k,N)+S2 {q, t;j,k,N) 

+ S3{q,t;j,k,N) , (4.5) 

where 

S,{q,t;j,k,N) = e-^^\f'-jit)\ (4.6) 

is the contribution ignoring j< , j> , thus ignoring all end 
effects. <S2 corrects <Si for the constraint j> > j< 

52(g,t;j,fc,7V) = [e(j>-j<)-l] e-9'l^-i(*)l , 

(4.7) 

and S3 takes the newly created parts of the tube into 
account 
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S3{q,t;j,k,N) = Q{k-j{t)) e(i<-fc - 1) e(j>-j<) 



,-q\2u-k-j{t))_.-q'{k-jit)) 



+ e(j-(t)-fc-i) e(j<-j(i)-i) e(j>-j<) 


■^-f{2u-k-m) - e-fim-k)' 


+ 0{k-j{t)) e(j(0-j> - 1) eo>-j<) 


e-9'(fc+j(i)-2j>) _ ^-f{k-j{t))- 





+ e{jit)-k-i) e{k-j>-i) e{j>-u) \e-f{k+jit)-2j>) _e-f{j{t)-ky . 



(4.8) 



These expressions are formally exact, but, as pointed 
out in the previous section, to cvahiatc S2 and S^, we 
have to construct an approximation for the simultane- 
ous distribution of n{j,t), j<, and j>. Si could be 
evaluated with the exactly known distribution of n{j, t) 
(Eq. (3.3)). However, being interested in the universal 
features of the model, which only show up for larger 
times, we use the Gaussian approximation (3.11) and 
ignore the discreteness of the elementary hopping pro- 
cess. We also will take the chain as continuoiis, in the 
evaluation replacing segment summations by integrals. 
A priori these simplifications might influence the short 
time behavior, but in practice they are found to have 
no measurable effects. For a check, we numerically have 
compared the continuous model to a fully discrete eval- 
uation. For the properly normalized coherent structure 
function Sc{q,t; M, N)/ Sc{q,0; M, N), the difference for 
all times, including the microscopic range, is found to be 
of the order 10~^ and thus negligible. 



B. The contribution Si{q,t;j,k,N) 

Combining equations (4.2) and (4.6) with the defini- 
tion (3.3) of Pi, we find 

+00 _„ 

Si{q,t;j,k,N)= J2 e-^\^-^-^s^\Vi{n;j,t). 



(4.9) 

With the Gaussian approximation (3.11) for Vi{n,j,t) 
and with n taken continuous, this expression is easily 
evaluated to yield 



Siiq,t;j,k,N) 



(4.10) 



e^^"^ erfc (Q + A) + e'^^^ erfc (Q - A) 



where 



Q = q^is ^/poAi{j,t) 
A k - j 

2£sVPoAi{j,t) 

Note that Eqs. (3.2) and (3.7) imply 

^PoAi{j,t)^({rj{t)-rjiO))^) , 



(4.11) 



so that in the result (4.10), q and k — j are measured 

relative to the mean displacement of segment j. Even 
though the Gaussian approximation from its derivation 
holds only for poAi 1, it for poAi 0, i.e., t — > 0, 
reproduces the exact static behavior of our model: 



Siiq,0;j,k,N) = e 



(4.12) 



As will be discussed in Sect. V.E, this is an important re- 
quirement for any theory of the dynamic scattering func- 
tions. 



C. Distribution function for S3{q,t; j, k, N) 

In view of the symmetry of the chain under reflection 
j ^ N — j, the last two terms in Eq. (4.8) for 53, when 
summed over j and k, yield contributions identical to the 
flrst two terms. We therefore can restrict the analysis to 
the first terms, which involve the distribution function 
referred to in Sect. III.B (Eq. (3.20)): 



For instance, in terms of this distribution function, 
the contribution to the coherent scattering function 
Sc{q, t, N) of the first term in Eq. (4.8) reads 



J2 E ®ik-j-lnj)Q{l,nm-k-l)V^ij{nm,nj;t) 

i,k nm,nj 

-q^{2isnm - k-j - IsUj) _ ^-f{k - j- IsTij) 

We now construct an approximate expression for V^^^ j , 
based on random walk theory. We first present the essen- 
tial steps of our approach and discuss the approximations 
involved thereafter. Some details of the calculations are 
given in Appendix A. 

We introduce the auxiliary variable no = — n(0, t) and 
write 

-t-00 



no=—oo 
-t-00 



= E 



,o,jinm,no,nj;t) . 



(4.13) 



7 



Recall that n(0, t) is the change in the occupation of 
reservoir within time t. We then factorize according 
to 

'Plnl,o.jinm,no,ny,t) (4.14) 

"^max,o('^m, no; t) 7^2 (»o, nf,OJ, t) 
Pi(no;0,i) 

IT'] 

Here V^J^ f^{njyi,no;t) is the simultaneous distribution 
of nmax(0, i) and —n{0,t), with the constraint j> > 
j< taken into account. V2ino,nj;0,j,t) is the simul- 
taneous distribution of —n(0,t) and n{j,t), so that 
V2{no,nj;0,j,t)/Vi{nf);0,t) is the conditional probabil- 
ity to find Uj = n{j,t), once tiq = —n(0,t) is given. A 
rigorous expression for V2 was given in 1^, Sect. III. Here 
we again use the Gaussian approximation ( , appendix 
C, Eq. (C.7)). Using also the Gaussian approximation 
(3.11) for Pi{no;0,t) we find 



'P2{no,nj;0,j,t) 
■ exp 

where 



2(1 -a2) 



27rn2(j,i)) ^'^ (I 



-1/2 



no 



(Ai(0,<) Ai(j, i))i/2 
with [cf. Ref. If, Eq. (A. 12)] 



(4.15) 

(4.16) 
(4.17) 



Pi 

N 



N 



1 

ON 



3 

2N 



2N 



1 



N-l 

2N ^ 

K=l 



cos(#)cos(^(j + i)) 



(4.18) 



y2N) 



A\ and are given in Eqs. (3.5), (3.6), respectively. 

This introduces the parameter function a — a(j,t). It 
measures the coupling of the motion of bead j to the 
motion of chain end 0. If it vanishes, the conditional 
probability (4.15) reduces to Vi{nj; j,t). This happens 
for <C (see Ref. 1^, Eq. (A. 16)). The maximal value 
of a is 1, which is approached for j — > and all i, or for 
t > T2 and aU j. In the latter hmit, Eq. (4.15) yields 
Zj — zq, and the motion of all segments is rigidly cou- 
pled to the motion of the end-segment. In this limit, we 
thus recover the basic assumption of the primitive chain 
model. 

To find an acceptable functional form for 
Pj^a^ Q(nm, no; t), we replace the correlated stochastic 
process n(0,s) by a random walk n'(s) on the integer 
numbers, with hopping rate p'. We have to consider 
walks that start at n'(0) ~ 0, end at n'{t) = no, and 



attain the maximal value rim > no for some s S [0, t]. To 
take care of the constraint 

i(j> -J<) = N'- n,„ax(A^, t) - n„,ax(0, t) > , 

(4.19) 

we restrict the walk n'(t) to the interval [rim — N' + 1, n^], 
where N' is the greatest integer less than N/£s, and we 
use absorbing boundary conditions. This amounts to the 
assumption that a particle entering the chain from the 
reservoir at chain end N is transfered immediately to the 
reservoir at chain end 0. This assumption is in the spirit 
of the primitive chain model. 

With these simplifications, V^J^ q can be calculated as 
sketched in Appendix A. Our result reads 



'PmLoinm,rio;t) 



e n 



-\-OC 

E 

U——00 



|no|-t-no 
2 

(i^ + l) 



e i^N'-rim- 
vN' 



|no|-no 
2 

Tm_ _ no 

/71 2V7< 



exp 



vN' 
no 



■ exp 



no 



no 



/j7t 2y/j7t 



(4.20) 



It is valid for p't ^ 1 and ' ^ 1 , which is the region of 
interest here. We now note that for p't ^ 1, the relations 



n'2(t) = 2p't 




max n' (s) 

se[o,t] 



hold, and we use these relations to eliminate the factors 



no 



= V2 



no 



= V2 



V2 J^^ = V2y 
V TT 



/jyt 



N 



IT ZMr^ 



= V2N 



In the last line, we used n^, since N' is introduced via 
the constraint (4.19). Furthermore we write the prefactor 
in Eq. (4.20) as 



/tt p't 



-1/2 
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and we again treat zq and y as continuous variables, 
which for p't ^ 1 is a vahd approximation consistent 
with our derivation. With these substitutions, Eq. (4.20) 
reads 



2 nmax(0,i) 



dno dUrn V^^^n 



^ f^max(0) 

N = N/c , 



2 y V 2 

I (ly+l) {2uN+ 2y-zo) exp 



+ 00 



i^——oo\. 



dzo dy 



{2iyN+2y-zoy 



with 



C = c(t) ^ ^j-ts nmax(0, t) 



(4.23) 



(4.24) 



(4.25) 



I' {2i'N+zo) exp 



{2i^N+zay' 



(4.22) 



For nniax(0,t), the approximation (3.16) is used. Up to 
the factor -^/f , the parameter c gives the distance up 



to which the tube has been destroyed on average. We 
now use expressions (4.22), (4.15), (4.14) to evaluate Eq. 
We use this result which has been derived for a random (4.13), where the sum over uq has to be replaced by 
walk, also for the correlated process n(0, s) by reinter- the integral over zq. Some exercise in Gaussian integrals 
prcting the variables, zq is given by Eq. (4.16), and yields our final result: 



+ 00 



dn,dn™7'£^(n„,n,;t) = ^e(2/)e(iV-y) 



(4.26) 



^mixj(y'2:,a) 



+ 1) exp 



(2auN+2ay-z)'^ 



V exp 



{2ai'N-zf 



with the notation 



a{2avN+2ay~'z) 



^ , az-y + a2y + a2vN\ I az-y + a2y + a2vN + N\ 
eric I — I — eric I — 



/ {az~y+a2y + a2vN)'^\ 
exp — exp 



/a2 / 
{az~y + a2y+a2vN+N)'^ \ 



02 



02 



a{2aiyN-z) 



erfc 



az — y + 



a2vN\ 



0.2 I 



erfc 



az-y + a2i'N+N 



/a2 



a2 = 2(l-a^) 



exp 

2\ 



(az-y+a2i'NY\ ( {az-y+a2i^N +Nf 
— - exp 



02 



02 



(4.27) 



Clearly, the approximations inherent in our construc- 
tion of ■P^xj need some justification. Steps hke the 
replacement of discrete by continuous variables are well 
justified, since we need the result only in a time and 
chain-length regime where a continuous chain model is 
valid. The problematic steps are the factorization (4.14) 
of 'PI'^I n , and the calculation of the functional form of 

max, U.J 

(T) 

'^max random walk theory. 

Technically, the factorization (4.14) serves to reduce 
the problem to the treatment of the single stochastic 
process n(0, s). It clearly is justified for large times, 
t ^ T2, where n{j,t) is firmly bound to n(0,t), and 

(T) (T) 

where V^J^ j becomes equivalent to Pmax,o- ^'^^ shorter 
times t ^T2, it assumes that —n(0,t) is a good measure 
of ?^max(0,^), which is certainly incorrect, in particular 
for i ^ T2. However, for t ^ T2, end effects influence 



only a small part of the chain. As can be seen from 
Eqs. (3.16), (3.7) and has been exphcitly worked out in 
Ref. g, E q. (5.28), ^w(M) for t < T2 behaves as 
nmax(0, t) ~ (pt)i/4 - N^/^{t/T2y/'^ < N. Since the co- 
herent scattering function sums over all segments, it for 
such times is dominated by the motion of interior seg- 
ments not influenced by end-effects and governed by the 
distribution function Vi{n;j,t). It is easily verified that 
in the appropriate limit N = N/c{t) 00, the distribu- 
tion function (4.26), when integrated over n,„, reduces to 
the Gaussian approximation for Vi (n; j, t). (Note that in 
this limit only a part of the 1/ = contribution to the 
sum in Eq. (4.26) survives.) 

We now turn to our construction of P^^j^^ 0(71^, no). 
Our treatment of the constraint (4.19) should be ade- 
quate, since this constraint is relevant only for times of 
the order of the reptation time, i « T3, where the inter- 
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nal degrees of freedom of the chain are irrelevant. Fur- 

(T) 

thermore, V^J^ q by construction obeys the constraint 
no < rim, and integrating over no, we find a d istribution 
with the desired first moment n„i — nniax(0, t). Integrat- 
ing over Um, we find the correct (Gaussian) distribution 
of TiQ. With these three important features guaranteed, 
we may hope that we have found a good approximation 
for the distribution function of the full correlated process. 

To summarize, our construction interpolates among 
two limits where the full dynamics reduces to that of 
a single stochastic variable. For t ^ T2, the motion of 
individual segments governed by Ti{j, t), is essential. For 
t^T2, the parameter function a tends to 1 and the in- 
ternal motion becomes irrelevant. Furthermore, the cor- 
relations of the stochastic process n(0, s) have died out. 
We thus are concerned with a single uncorrelated pro- 
cess n(0,s), as in the primitive chain model. Smoothly 
interpolating among these limits, we may hope to have 
found a good approximation also in the crossover region 
t « T2. Indeed, as shown in Figs. 10 and 11 of Ref. 
a simplified version of our theory almost quantitatively 
fits with Monte Carlo data for the motion of individual 
segments. Furthermore, as will be illustrated in Sect. 
VIII and in more detail in Ref. |lj] , our theory quantita- 
tively accounts for data for the coherent scattering func- 
tion Sc{q,t; M, N). The agreement is equally good for 
the total chain (M = N) where tube length fluctuations 
are very important, and for an interior piece (M < N) 
where tube length fluctuations are irrelevant. 



D. Distribution function for S2{q,t; j, k, N) 



To evaluate S2{q,t;j,k,N) (Eq. (4.7)), we need the 
distribution function 



rPinjit) = [e(j> - j<) - 1] ,5„^,„(,-,) 

00 

= E ^-aL- i^rn , n, ■,t)~V, (n, ; J, t) . (4.28) 



The first part can be determined by integrating T'j 
(Eq. (4.26)) over n^. Eqs. (4.26)-(4.28) thus yield 



(T) 
max j' 



, +00 



J2 VY\z,,a) (4.29) 



U——00 



[zj ,a)=^ dy v'^^l^j {y, z^ , a) - 6^,0 exp(-zJ/2) . 
"'0 

(4.30) 

The 2/-integral can be evaluated analytically and we quote 
the result, which is useful for the numerical evaluation of 
^2, in Appendix B. 



E. Result for the tube conserving contribution to 
the coherent structure function Sc^^ [q, t; M, N) 

We consider the scattering from the M + 1 central 
beads (Eq. (3.18)) and write 

SP {q, t- M, N) = Si {q, t; M, N) + S2 (<Z, t; Af , N) 

+S3(q,t;M,N) , (4.31) 

where the Si{q,t; M, N), i ~ 1,2,3, are the contribu- 
tions Si{q,t] j,k, N) (Eq. (4.5)), summed over j and k. 
The superscript '^^^ again recalls the constraint (4.1): 
J> - J< > 0. 

Due to this constraint the relation 



S,{q,t-M,N) = SP{q,t;M,N) 



(4.32) 



in general holds only for t ^ T3. However, for large q 
such that (f'R^ ^ 1, contributions in which the tube has 
been destroyed, contribute neghgibly to Sc{q, t; M, N), so 
that Eq. (4.32) in this limit holds for all times. 
Consider now the first contribution. 



Si{q,t;M,N) = 



dj dk Si{q,t;j,k,N) 



Using Eqs. (4.10) and (4.11), we can carry out the inte- 
gral over k to find 



Si{q,t;M,N) = ^ 



dj 



e^AiQ+Q^ erfc(Q-)-Ai) - e-^^''^+'^ erfc(Q-Ai) 
„ g2A.Q+Q^ erfc(Q + A2) + e-2A2Q+Q^ erfc(Q- A2) 



2 erfcA2 - 2 erfcAi 



(4.33) 



where 



Q^fisVpoAiij,t) 
1 /N + M 

1 fN-M 
2L 



Ai 
A2 



,j] ipoA,{j,t)r^/^ . (4.34) 



The remaining integration in general must be done nu- 
merically, due to the j-dependence of Ai{j,t). 

The integral over k can be carried out also in 52,53. 
We introduce the notation 



q^fc; j=j/c; M^M/c 
1/2 



5=^(n2(^-f j,t)) 
c 



(4.35) 
(4.36) 



where the parameter c has been defined in Eq. (4.25). 
With due regard of the 8-functions, a straightforward 
calculation yields 
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S2{q.t-M,N) 



Q Jo 



M 



s: 



(2) 



(4.37) 



(2) 



/2n 



E 



M 

dz 



2 - 



-g(M-^) 



, a 



(4.38) 
(4.39) 



o 2 i-M 

53 (g, t ; Af, TV) = ^ / (si'^+ 5f ) + 5f ) + S: 
Q Jo ^ 



(4.40) 



^(N+M)/2 i-l 



5. 



(2) 



s. 



(3) 



dyy I dz 





2 1 



1 



E^. 



Jo Jo V^TT 



max, J 1 o 



/27r 

N - M z+j 



^ l.(N-M)/2 j-l 

dyy dz 
Jo 



max, J 



N + M y M-j ^ 



(4) 1 " ' 



(4.42) 



(4.43) 
(4.44) 



The prefactor of 2 in Eq. (4.40) accounts for the last 
two contributions in Eq. (4.8). We note that these resuhs 
depend on time via the parameters a, c, and b. From its 
definition (4.36), the parameter function b — b{j,t) mea- 
sures the motion of an arbitrary segment relative to the 
motion of the end segment. It is weakly dependent on j 
and tends to 1 for t ^ T2. 

In ^2 and ^3, one more integration could be done an- 
alytically which, however, only blows up the number of 
terms without leading to any simplification. Due to the 
dependence on the segment index j implicit in a = a(j, t) 
and b — b{j^t), an analytical evaluation of all integrals 
is possible only in the limit t ^ T2 where a 1 and 
b —f 1. In general, we have to resort to numerical evalua- 
tion. In this context, we may note that the summations 
over V for i ^ T3 converge rapidly, so that in the range 

where sP{t ^ 0)/5c(0) exceeds lO'^, we never need to 
go beyond < 4. 



V. DYNAMICS WITHIN THE INITIAL TUBE 

In a time region where end effects are unimportant, the 
results of the previous section can be simplified. In pre- 
cise terms, the neglect of end effects amounts to consider- 



ing a subchain of length Af , in the center of an infinitely 
long chain. We here concentrate on this particular limit 
and compare our results to those derived for a Rouse 
chain in a coiled tube. 



A. Results of the reptation model 



In the limit N — > 00, with j = j — and k ^ k — ^ 
fixed, only the contribution Si to S'^'^\q,t; j, k, N) (Eq. 
(4.5)) survives. Furthermore Ai{j,t) (3.5) simplifies to 
[see Ref. [|, Eq. (4.11)] 



TT 



dz 

7^ 



At 



(5.1) 



independent of j. (Recall the definition t — pt.) 
Si{q,t;j,k,N) (Eq. (4.10)) takes the form 



Si{q,t;j,k) 



(5.2) 



e^^^ erfc (Q + A) + e^^Ag ^^^^ (g 



where now 
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k-j ft 



1/4 



-1/4 



(5.3) 



Integrating over j and fc, we find for the normalized dy- 
namic structure function 



Sciq,t;M, oo) 



Sc{q,t;M, oo) 
AP D(q^M) 



l-e 




(5.4) 



- 2 erfc(Q) 



We recall that D{x) = ^(e ^ — 1 + x) is the Debye 
function, and note that (fAI can also be written as 



(5.5) 



where Rg{M) is the radius of gyration of the subchain. 
In Eq. (5.4), we introduced a new time scale 



^4 = ^(^1po)-^M^ 



which is of the order of the time the subchain needs to 
leave its original part of the tube. This interpretation 
is obvious from the results on segment motion quoted at 

the end of Sect. II: (^{rjifi) - rj{0))^'j ~ fl^"^ - M for 

T4 ^ T2. Also the variable Q can be expressed in terms 
of a time scale: 



Ur 



1/4 



(5.7) 



then find the following limiting behavior in the various 
time regimes 



M2 



(5.9) 



il-elllfi^ 

' (<?2_R2(Af))2 

' (A) 



(*) 

: \ 1/4 



1/2 



for t < Ta 



for < t < T4 



Sciq,t;M, 00) 
M2 



6 Tj 
V^q^UM) [ i 



1/4 



for r4 < < . 
(5.10) 



Note that our results depend only on macroscopic pa- 
rameters q^R^g{M) and i/Ti^ which absorb any reference 
to the microscopic structure. 



B. Comparison with de Gennes' results for Rouse 
motion in a tube 



In Sect. II. B, we recalled de Gennes' results |T^, de- 
rived for one dimensional Rouse-type motion in a coiled 
tube. In the derivation, the relation q^R^g » 1 S> q^P 
was assumed, with ^ being the average segment size of 
a Gaussian chain. Of interest here is the 'local' term 
S^^\q^t) (Eq. 2.12). A glance to the derivation shows 
that it implicitly exploits the limit considered here: a 
subchain (of length M) in an infinitely long chain. Com- 
bining Eqs. (2.10) to (2.14) we thus find 



(5-6) 5rfG(g,t;M,oo) 
M2 D{q^M) 



= 1 - r 



9 erfc 



T, 



const 



t/Tg 

(5.11 



(5.12) 



Clearly this expression differs strongly from our result 
(5.4). It leads to very different asymptotics: 



1674 



iq'RUM)y 



(5.8) 



It needs the time Tq before the distance diffused by the 
subchain, can be resolved by scattering of wave vector q. 
For comparison with previous work, we concentrate on 
q^Rl{M) > 1, so that fg < fi and D{x) « 2/x. We 



SdG{q,t;M, 00 



D[q^Rim 

1/2 



2iVe 

M 



^ (£) for t « Tg 

1/2 



for t > Ta 



(5.13) 



Furthermore the scaling with q^ ^ M, and t is quite dif- 
ferent. 
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C. Closer inspection of Rouse motion in a tube 

The derivation of Eq. (5.11) in Ref. ||l^ involves some 
approximations, which greatly simplify the analysis but 
are not really necessary. For a general test of the valid- 
ity of the model, we therefore have repeated the analysis 
without these simplifications. The analysis is sketched in 
Appendix C. For the further discussion, we here quote 
the result for the scattering from a given pair of beads 
(Eqs. (C.16) - (C.18)). 



SRT{q,t;j,k) 



|e2^^ erfc(g + A) 



(5.14) 
e-2AQ erfc(Q- A)} 

1/2 



A = 



6 \3-~k\ 



\J - k\ 



: erfc , 



|J - fc| 



-1/2 



(5.15) 



(5.16) 



Here 70 is the segment mobility of the one-dimensional 
Rouse model. 

Clearly the structure of Srt (Eq- (5.14)) is identical 
to that of our result Si (Eq. (5.2)). The difference is in 
the quantities Q, A (Eq. (5.15)), compared to Q, A (Eq. 
(5.3)). We note, however, that the relation 



2Q/^^qH'\j -k 



q^^^\j-k\ = 2Qk 



holds. Recall that the mean segment size ^ of a Gaus- 
sian chain, which is asymptotically equivalent to a chain 
with fixed segment length £0, obeys = £9/6 (in three 
dimensions). 

To analyze the difference among the two models, we 
first consider the static limit, t = Q. Si reduces to (recall 
the definition (f = = q^t"^) 



Si{q,0]j,k) = e 



(5.17) 



which is the exact result. The result for But can be 
written as 



Even if we ignore the terms q^i'^N^/Q ^ 1, taking them 
to be irrelevant micro-structure effects, this result does 
not reduce to the exact form (5.17). To recover this form, 
we rather consistently have to take the limit q^l^ — ^ 
with q'^£'^{j — k) fixed, i.e., \j — k\ = \j — k\ ^00. 
This just demonstrates that a one-dimensional Gaussian 
chain, folded into the three dimensional random walk 
configuration of the tube, does not yield the exact distri- 
bution of a three dimensional Gaussian chain. In other 
words, the model of a Rouse chain in a tube violates the 
equilibrium initial conditions by micro-structure terms 
on scale \j — k\ ^ N^. 

For the dynamics, this discussion implies that the 
model gives a wrong estimate for the number of wiggles 
in the initial configuration. To eliminate this effect of 
unphysical initial conditions, we have to take the same 
limit q'^i^ with q^(.'^{j — k) fixed, also in the full time 
dependent expression (5.14). To facilitate the discussion, 
we rewrite Eqs. (5.15) in a form which exhibits the fixed 
combination of variables q^(-'^\j — k\. 



Q^qi 




'■^ I y \/ £2 y I ^2^2^2f^ 



1/2 



—J—qH'^lj -k\ 



q^e\j^k\ + q'e\i2:§tg 



qH^\j~k\ 



(5.19) 

-1/2 



The function g{z) obeys the relations 



5(0) = ^ 



Now the limiting result for Srt sensitively depends on 
the way we scale the time. We first consider times such 
that q^y/t stays finite upon taking the limit 0: 



qH\l2^t^q-^, a>0 



We then find 



SRT{q,0;j,k) 



exp 



-q^e\j-k\[l-q^e^^ 



1 



1 



exp 



2 

•erfc 



erfc 



N, 



■\j~~k\ 



q r\j-k\ i + q^i 



— \j-M 



,2 «2 



2 ' 
6 
6 



(5.18) 



g ^ 0, A -> 00 

and recover the static limit (5.17). Indeed, for such times 
the scattering cannot resolve the internal motion: 



q' ({V,it)-V,m)^q't 



Effects of internal dynamics can be seen only for times 
such that q'^Vt diverges: 
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a > 



Then the contribution proportional to g dominates the 
square brackets in Eqs. (5.19), and furthermore the ar- 
gument of g tends to zero. Thus 



1/4 



-1/4 



(5.20) 



which is the same functional dependence on t and |j — A;| 
as that of Q or A (Eq. (5.3)). The limiting expression for 
Srt{q, t; j, k) becomes identical to our result Si{q, t; j, k) 
if we identify 



2*2 



18 £2 



£2 



t 



;po) t 



(5.21) 



(5.22) 



In summary, we have found that the model of a Rouse 
chain in a tube is equivalent to the reptation model only 
in the hmit q'^i'^ with q^i^{j — k) fixed. Outside this 
limit, it exhibits an unphysical relaxation of nonequilib- 
rium initial conditions. 



D. Relation among the microscopic parameters of 
the different models 

As a byproduct of our analysis, we with Eq. (5.22) 
have found a relation among the microscopic parameters 
of our model and those used in more standard Rouse type 
modeling of chain dynamics. Analyzing in Sect. VII the 
relation of our model with the primitive chain model, we 
will find as an additional result: 



Ne To 
6 P 



Combining Eqs. (5.22) and (5.23), we find 



70 

2^2 



t = t 



(5.23) 

(5.24) 
(5.25) 



We now can give a quantitative definition of the equili- 
bration time T2, which we identify with the Rouse time 
of a free chain of N segments: 



To 



2 /2 

4(^+1) - 

TT^ 70 



With Eq. (5.25), we find 



T2=pT2 



(5.26) 



E. Implications for the coherent structure function 

The artifact of the model of a Rouse chain in a tube 
concerns only small parts of the chain of the order of the 
tube diameter and thus might be thought to be negligi- 
ble. In the static structure function, the error sums up 
to a term of order N, small compared to Sc{q,0;N) — 
N^D{<f'N). Since, however, for t ^ T2 the time depen- 
dence of Sc is weak, even such a small effect is relevant. 
Indeed, it can greatly change the picture. To illustrate 
this, we in Fig. 3 compare our result (5.4) for the normal- 
ized coherent structure function Sc{q, t; M, 00) to the re- 
sult found by integrating SRT{q,t;j,k) (Eq. (5.14)) over 
<j < ^ a.nd <k < ^. To relate the models, 
we used relations (5.21), (5.24), and (5.25). To include 
also de Gennes' approximate form (5.11), we used the 
large value q^R^{M) = 50. We note that de Gennes en- 
forced the correct i = behavior by artificially subtract- 
ing his result for Sjit- We furthermore note that repeat- 
ing his calculation in our notation, we found t/Tq = qH. 
The remaining parameter (."Ipo = 1.23 has been taken 
from our previous work |^ . 

Fig. 3 shows that the effect of the artificial initial con- 
ditions can be quite large and dies out only slowly. It 
extends up to the Rouse time of the subchain considered. 
This result is generic. For longer chains, the amplitude 
of the effect decreases for the normalized structure func- 
tion, as expected for a micro-structure effect, but the 
range stays of order T2(M). This is obvious since a non- 
equilibrium initial condition generically will relax only on 
time scale T2. As a side issue, we note that de Gennes' 
approximation (5.11) agrees quite well with the shifted 
form of Srt- 

To close this section, a general remark on micro- 
structure corrections for the dynamics may be appropri- 
ate. Our result shows no such corrections, which would 
give rise to an additional 1/Af -dependence in Eq. (5.4), 
which is not in the form of the scaling variables (fM and 
t/T4. Thus our model succeeded in singling out the uni- 
versal aspects of reptation dynamics. This, however, does 
not imply that (non-universal) terms yielding some ad- 
ditional 1/M dependence could not show up for a micro- 
scopically realistic model, which takes the details of the 
microscopic motion into account. But we want to stress 
that any model first of all has to yield the correct static 
structure function. Otherwise some unphysical relax- 
ation will infiuence the dynamics. Such results can safely 
be trusted only in a range where Sc{q, 0; N) — Sc{q, t; N) 
exceeds the error in Sc{q, 0; N). 

We finally note that here we have been concerned ex- 
clusively with the reptation aspect of the dynamics, mod- 
eled as one dimensional Rouse motion in a tube. This is 
to be clearly distinguished from three dimensional Rouse 
motion among fixed entanglement points, as treated by 
Des Cloizeaux [|l^, for instance. The latter model is 
concerned with motion in melts for 'microscopic' times: 
t<To. 
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VI. ANALYSIS OF COMPLETE TUBE 
DESTRUCTION 

In this section, we derive an integral equation extend- 
ing Sc{q,t; M, N) to arbitrarily large times (subsection 
A). Basically it is an equation for S{q,t;j,k,N), which 
incorporates S^'^^ {q, t; j, 0, N) as inhomogeneity. To cal- 
culate Sc, we need to sum this inhomogeneity over j. We 
construct this function in subsection B, following the ap- 
proach of Section IV. The kernel of the integral equation 
involves some distribution function which is calculated 
in Subsect. C with the help of the mean hopping rate 
approximation. Quantities like the probability density of 
tube destruction at time t, which can be derived from 
this distribution function, are discussed in Subsection D. 
The numerical evaluation of our results for Sc{q, t; M, N) 
is deferred to Sect. VIII, after we have shown that our 
theory in the appropriate limit yields the results of the 
primitive chain model. 



A. Derivation of an integral equation for the 
structure function 

Up to now, we only considered stochastic processes for 

which some part of the initial tube still exists at time t. 
To get rid of this constraint, we have to deal with situa- 
tions as shown in Fig. 4: at time to, < to < t, the chain 
leaves the original tube, which means that the remainder 
of the original tube is the single point rjg{0). This point 
is occupied by a chain end. The rest of the chain has 
found a completely new configuration. 

Let V*{jo,to I 0) or P*{jo,to \ N) be the probability 
that the tube is finally destroyed at time to, the last point 
rj(,(0) being occupied by chain end or N, respectively. 
We assume that V* does not depend on the initial config- 
uration, which should be satisfied except for rare extreme 
cases. 

We then can write the full time dependent scattering 
(3.17) from a pair of beads j, k as 

S{q,t;j,k,N) = S^^\q,t;j,k,N) 

t-l N 

^Y^H [^*(io,to I 0) S{q,t,j,k\ io,io,0) 



m = 0,7V 



(6.2) 



to=l jo=0 



+V*{jo,to I N) S{q,t,j,k\ jo,to,N) 



where S{q,t,j,k \ jo,to,m) with m = O.N denotes the 
scattering with tube destruction specified by jo, to and 
m. We now factorize S{q, t,j, k \ jo, to, m) according to 

S{q,t,j,k I jo,to,m) 



„iq(r,(t)-rfc(0)) 



ja-ta,rn 



M{vj{t) - Vm{to)) Jq(rjo(0) - rfe(O)) 



Jq(rj (t) - r„(to)) 



where the second factor in the last line is a purely static 
average. We have exploited Vm^o) = ^jo(O)- This factor- 
ization should be well justified, since the chain at time 
to has attained a completely new internal configuration. 
Now the first factor in the last line of Eq. (6.2) equals 
5(g, t—to;j, m, N), whereas the second factor is the static 
structure function cxp[— \ jo — k \]. Combining Eqs. 
(6.1), (6.2) we thus find 

S{q,t;j,k,N)=S^^\q,t;j,k,N) 

t-i 2 
+ E E [P*ijo,to I 0) e-9 lio - k\ s{q,t-to;j,0,N) 

to=l 30 

+V*{jo,to\N) e-9'bo-/c| S{q,t-to;j,N,N) 
Reflection symmetry along the chain implies 

r*{jo,to\N)=V*{N-jo,to\0) 
S{q,t;j,N,N)=S{q,t;N-j,0,N) , 



so that our result takes the form 

S{q,t:j,k,N)=S'^^\q,t:j,k,N) 



(6.3) 



-EE^*(jo,to|0) 

*o=l jo 



-9'\jo-k\ s{q,t-to;j,0,N) 



-f\N-jo-k\ 



S{q,t-to;N -j,0,N) 



We now sum j and k over the central piece of the chain 

to find 

S,iq,t;M,N)^ SP{q,t:M,N) 

t-l ( N (N+M)/2 

+2 E E E ^*(io,io|0)e-9'lio-fc| 

«o = l \3o=Ok=(N-M)/2 

■SE{q,t-to;M,N) , (6.4) 



(6.1) where 



(N+M)/2 

SE{q,t;M,N)= ^ S{q,t; j,0,n) . (6.5) 

j=(iV-M)/2 

For Se, Eq. (6.3) yields 
SE{q,t;M,N)=s'-P{q,t;M,N) 

t-l ( N 

+ E E T*Uo, to I 0) (e-«'io + e-fiN - jo)) 

■SE{q,t-to;M,N) . (6.6) 

With time (and segment index) taken continuous, this is 
the basic integral equation of our approach. We note that 
it is of Volterra-type and therefore has a unique solution. 
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B. Expression for S)^ ' (q,t; M, N) 

The tube conserving contribution to SE{q,t;M,N) is easily found from the results of Section IV. Following Eq. 
(4.31), we write 



SE,i{<l:t;M,N) 



{q, t; M, N) = 5^,1 (g, t; M, N) + SE,2{q, t; M, N) + SeMq, M, N) . 
djSi{q,t;j,0,N) 



N-M 



r(N+M)/1 i-oo ^-qz I 

SE,2{q,t;M,N) =c I dj dz—^-j=Y, 



SE,z{<l^t-M,N)=c [ 



{N-M)/2 
(JV+M)/2 ^ ^ 



(6.7) 

(6.8) 
(6.9) 





/•OO 








e 









-q{2y + z) _ -qz 



(6-10) 



+ dzy 



-qy{2 -z) _ p-qyz] _J_ 



yz-j 



2^ ^ max,, . y> ^ ,a 



The notation is the same as in Sect. IV (see, in particular, Eqs. (4.27), (4.35), (4.36)). 

C. Expression for V*{jo,to | 0) 

To construct an expression for V*{jo,to | 0), we again use random walk theory, closely following the derivation of 
V^J^ '^^ Sect. IV. C. The calculation is sketched in Appendix A. It yields the result 



dtoV''ijo,to\0) = j^inp'{to-l))-'^^ ^ 



+ ■00 



v——<yo 



2 (z.7V + jo/2)^ 



(io - 1) 



exp 



(i/iV + jo/2)^ 
■ (to - 1) 



(6.11) 



We again express to — 1 ~ in terms of the maximal excursion Um (Eq. (4.21)) and identify rim with nmax(0,to). 
This yields the replacement 



4Vp'(t0 - 1) ^ ^ 4nma^(0,to) = 



resulting in 



N 



, =^x/2 7V = y2- 

, Vp'(*o - 1) c 

dig ^ 2 — 
to - 1 c 



Jo 



4Vp'(to-l) 

4 V^P'(t0 - 1) ' V TT 



x/2 3o = ^^ 



- djo . 



With these replacements, we find 



2 dc 



+ 00 



djodtoV*{jo,to\0)^2J djo V 1 - 4{iyN + jo/2f exp -2(!.7V + jo/2) 

v TT C ■^-^ 



(6.12) 



I/— — OO 



To construct the kernels for the integral equations (6.4) and (6.6), we basically need 

dto r djo V^jo, to I 0) = 2 - V*{q, X, N) , 

Jo c 

where X = X/c, and where P*{q, X, N) is given by 



(6.13) 
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5Z ^ y "^Jo l-4[iyN + jo/2j cxp -2[i^N + jo/2 



+00 



J2 i^{fe^-^''+^/^ 



erfc 



2vN + X + q 



erfc 



2vN + q 



-J- 2vN + W - (2zyiV + 1-9) exp 

V TT V TT 



(2i/7v + i:)2 



In terms of V*{q, X, N), the kernel of Eq. (6.6) takes the form 



2 ICsiq, N) = dto djo V*Uo, h \ 0) 



e 9' JO _^g-r(w-io) 



2^ 

c 



TV, TV) + TV, TV) 



The kernel of Eq. (6.4) can be written as 



dc - . AN+M)/2 . 

2 — /Ce(9, TV, Tlf) = dto / (|?'o / dk V*ijo, h \ 0) e'l \Jo ' 

Q Jo J{N-M)/2 



= 2 



dc 



2V*{0,^,N)-2V*{Q,^,N) 

_ g9(JV-M)/2 j^-p*(^^ ^) _ TV)) + e-«(*-*)/2 V*{-q, TV) 

(P* (g, TV, TV) - V*{q, TV)) - e-«(*+*)/2 



„9(W+M)/2 



(6.14) 



(6.15) 



(6.16) 



D. Discussion of the probability density of tube 
destruction 



Comparing the present results with those of Sect. IV. C, 
we can verify the internal consistency of our random walk 
approximation. Yrom V^J^^ ^(nrmnj^t) (Eq. (3.20)), we 
can derive the probability tiiat a part of the initial tube 
still exists at time t: 

-00 N 
nj=—oo nm=Q 

A straightforward calculation starting from Eqs. (4.26), 
(4.27) yields 

p(^)(t) = l + 4E(-l)"i^erfc i^N 

= P(^) . (6.18) 

The probability that the tube is destroyed within time 
interval dto, corresponding to 

dc = c{to + dto) - c{to) , 



= _4^,/lTVy(-l)V e-i"'' 



(6.19) 



On the other hand, we can calculate this probability also 
as 



N 



2 dto E ^*(io,to I 0) =4 - P*(0,TV,TV) , (6.20) 

30=0 



where the factor of 2 takes the two chain ends into ac- 
count. It is easily verified that these two expressions are 
identical. Thus the following relation holds 



d_ 
dt 



N 



p(^)(i) = 2E^*(io,to|0). (6.21) 

30=0 



This identity guarantees the validity of the normalization 

M-N 
2 

SE{q = Q,t;M,N)= E ^ = M+l (6.22) 



can be calculated as —dc ^ V^'^'>. 



for all times. Prom the definition of , we have 
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j_N~-M_ 



(6.23) 



Substituting Eqs. (6.21) and (6.23) into the integral equa- 
tion (6.6), we find 

SEiO,t;M,N) = {M + l)V^'^Ht) (6.24) 

' dto (^^'^'(^o)) SEiO,t-to;M,N) . 

Partial integration together with 

SEiO,0;M,N) ^ M + 1 , p(^)(0) = 1 (6.25) 
yields 

0=/ dtoV'^^Hto) ^ SE{0,t-to;M,N) , (6.26) 
Jo c'f^o 



with only the trivial solution 



d 



— SEiO,t~to;M,N) = 
dto 



(6.27) 



Together with Eq. (6.25), this proves Eq. (6.22). The cor- 
responding analysis can be applied to Eq. (6.4), yielding 
the correct normalization 



SciO,t;M,N)= {M + lf 



(6.28) 



To get an impression of the time dependence of 
complete tube destruction, we in Fig. 5 a show 
N'^r*{0,N,N) ~ -f^V^'^){t) (cf. Eqs. (6.20), (6.21)) as 

a function of = 2{c/Nf = 7r(4 nmax(0, i)/iV)^ 

This choice of the variable is motivated by the rela- 
tion iV~^ ~ t/Tz, cf. Eq. (7.3). As we see, noticeable 
tube destruction starts at 2/iV^ w 0.1 and is essentially 
completed at k. 3.5. The variation of N'^V* as 

shown here, dominates the time dependence of the ker- 
nels (6.15), (6.16). It allows us to solve the integral equa- 
tion (6.6) for finite time i by a finite number of iterations, 
the result being exact within the numerical accuracy of 
our calculation. 

To close this section, we evaluate the probability that 
the initial tube finally is destroyed at the position of seg- 
ment Jo, with chain end being the last part residing 
in the initial tube (Eq. (6.12)). Fig. 5 b shows the de- 
pendence on jo/N for several values of N . As expected, 
for shorter times 2/7V^ < 1, chain end leaves the tube 
close to the other end {jo/N w 1). With increasing time 
the most probable point of final destruction slowly shifts 
to the center of the tube, but for times where the rate of 
the tube destruction is maximal, (corresponding to the 
maximum in Fig. 5 a), the shape of 'P*{jQ,to \ 0) is still 
quite un-symmetric. 



Obviously, the distribution functions considered here 
are closely related to the r.h.s. of Eq. (2.9), which is de- 
termined by the part of the original tube that is still 
occupied at time t (see Refs. fllJl^). In this context, it is 
interesting to note that Des Cloizeaux jl^ modified the 
expression (2.9) by replacing r = t/rd in the exponent 
by some more complicated time dependence, meant to 
take the local motion near an entanglement point into 
account. This modification is quite similar to our in- 
troduction of the quantity iV^ replacing t/rd. We note, 
however, that TV = N/c via the crossover behavior of 
c — c{t, N) takes end effects like tube length fluctuations 
into account rather than internal motion. 



VII. THE LIMIT OF LARGE TIME AND THE 
PRIMITIVE CHAIN MODEL 



A. Special cases 

If the time is large compared to the equilibration time 
T2 of the chain, our results simplify since the parameters 
a and b can be replaced by their limiting values 



1 



for t > T2 



(7.1) 



This implies that all segments experience the same curvi- 
linear shift, which is the basic assumption of the primitive 
chain model. Furthermore Ai{j,t) t/N for t 3> 72 (cf. 
Eq. (3.5)), and the parameter c (Eq. (4.25)) reduces to 



C = W - £snmax(0, t) 



(7.2) 



^ s— 1 ^ ^ 



Thus 



N \2T- 



1/2 



(7.3) 



becomes a direct measure of i/T^, where for brevity we 
introduced 



(7.4) 



as measure of the reptation time. With relations (7.1), 
all the integrals in Eqs. (4.37) - (4.44) can be evaluated 
analytically, resulting in a fairly lengthy expression for 
Sc {q,t; AI, N) as a sum of terms involving error func- 
tions and Gaussians. We here quote the result in those 

(T) 

limits, in which Sc becomes identical to the full scat- 
tering function Sc, which is the case for either short time: 
t/Ts ^ 1 or large wave vectors: q^Rg ^ 1. 
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1. Limit t/T'i with fixed Q — q^R^ 



We find 
Sc{q,t;N,N) 



t 

2T3' 



3/2 



(7.5) 



Recall that D{Q) denotes the Debye function. Of course, 
this limit can be attained only for an extremely long 
chain, since relation (7.1) implies T2/T3 0, i.e., N 
00. The result (7.5) shows that for such a chain relaxation 
becomes observable only for t ^ T2. Furthermore, with 
increasing Q, the time variation of Sc becomes rapidly 
insensitive to the scattering vector. 



2. Limit Q = q^R^ ^ 00 with fixed q^ and t/T-j, 

In this limit, our result reads 
Sc{q,t;N,N) 



00 

+4^(-l)^ 



- 1 r 



(7.6) 



' erfc 



71 



2 1 



which is the Poisson transform of 
S,{q,t;N,N) 



^ f;(2p-l)-2exp 



(7.7) 



-{2p-l) 



We thus recover the resuh of Refs. |2|jl3|, Eq. (2.9), pro- 
vided we identify 



t 



4 f. 



leading to 



{7- 



From Ref. j^], Eq. (6.19), we can take the relation of 
to the parameters of the underlying Rouse model, which 
in our notation reads 



Td 



(Replacement C, I/70, b ^ £o, 
Ref. §, Eq. (6.19).) Thus 



llN,, ksT = 1 in 



and Eq. (5.23) results. 



B. Proof of asymptotic equivalence to the primitive 
chain model 



Having recovered the results of Doi and Edwards for 
we clearly may ask whether for 
t ^ T2, the two approaches yield identical results irre- 
spective of q^R^. This is not obvious since formally the 



approaches are quite different. Doi and Edwards [[I3y2|] 
start from a diffusion equation for S{q,t; k, N). With 
the relation among model parameters established in Sect. 
V.D, this equation takes the form 



TV 9 92 



Siq,t;j,k,N)^0 



[pllpodt d,f_ 
This is amended by the initial condition 

Siq,0;j,k,N)^cM~f\j~k\) 
and boundary conditions 

lim ^ S{q, t; j, k, N) = q^Siq, t; 0, k, N) 
l\m ^ S{q, t; j, k, N) = -fS{q, t; N, k, N) 



(7.9) 



(7.10) 



(7.11) 



On the other hand, according to our theory, 
S{q^t; j, k, N) obeys Eq. (6.3), written in the continuous 
chain model as 

S{q,t;j,k,N) = S^^Hq,t;j,k,N) (7.12) 
+ j^dto djoV*ijo,to I 0) |e-«'l^«-'=l Siq,t-to;j,0,N) 
+ g-f\N-jo-k\ s{q,t-to;N - j,0,N)^ . 

The inhomogeneity takes the form 

S^^\q,t;j,k,N) (7.13) 
= dy I dz y2 V,y{y,z) F{y,z;j,k) 

Jo Jy~N ^ 

where 

'PAv^z) = J- lini Vi'^l (-'-'«) 

^ \ [(^ + 1) (2^7V + 2y-z) e-lt^i^'^N+^V'^f 



V (2vN — z) e ^ ' 



(7.14) 



The rescaling of the variables y and z serves to isolate 
the time dependence which now is contained in c only 
(cf. Eq. (7.2)). The function F collects all contributions 
contained in S\, S2, S3 (Eq. (4.5)), and is found to be 
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F{y,z;j,k) = e-'^"\''-^-'\ 

+e{k-j-z) Q{y-k) 



(7.15) 



^-f(2y-z-k-j) _ ^-f{k-j-z) 



+e(j + z-fc) Q(y-j-z) 



+e(fc ~ j + z) e{j - z + y- N) 



^-fik+j-z+2v-2N) _ ^-q^{k-j+z) 



+e{j-k-z) e{k-N + y) 



^-r{k+j-z+2v-2N) _ ^-q\j-k-z) 



.|e-9^bo-fc|-9^- ^ ^-f\N-io-^-f{N-^)^^ 

/■t pN 

+ dto djoV*{jo,to I 0) 
Jo Jo 

.|e-9^bo-fe| VS{q,t-to;j,0,N) 

+ e-f\N-jo-k\ vS{q,t-to;N-j,0,N)} . 

This is an integral equation for "DS which has only the 
trivial solution 



VS{q,t;j,k,N)=0 , 



(7.18) 



Here the last two contributions arise from the last two provided that the inhomogeneity vanishes. We first con- 



terms in Eq. (4.8), which a priori involve the distribution 
function 



sider the contribution VS^"^^ and note that in view of Eq. 
(7.2), D can be written as 



Interchange of the chain ends transforms this distribu- 
tion to 'P,^,^Lj('^m, -nj]t) (Eq. (3.20)) and imphes that resulting h 
we have to take j> = N — y and j{t) = j — z in the 
corresponding contributions to F{y, z; j, k). jq 

By construction, our form of S{q, t; j, k, N) obeys the "DS^^^ = dy 
initial condition (7.10). To derive Eq. (7.9), we apply the ^ Jo Jy 

operator 



X,--—- — 
c dc dp 



[ - ^My,z)]F{y,z;3,k) 



-N 



V ■ 



N d__cP_ 
P^odi'W 



(7.16) 



c dc 



- 'Pu{y,z) ■^F{y,z;j,k) \ . 



df 



(7.19) 



to the integral equation (7.12) to find 
VS{q,t;j,k,N)=VS^^\q,t;j,k,N) 



+ 



N 

pHpo Jo 



N 



djo'P*{jo,t I 0) 



rv 

VS^T)^^- dy dzFiy,z;j, 

^ Jo Jy-N 



F{y, z;j, k) (Eq. (7.15)) is a sum of terms which depend 
(7.17) on j and z exclusively via the combinations j + z ov j — z, 
respectively. Thus ^ is equivalent to and partial 
integration yields 



d_ 

dz 



F{y,z;j, k)-V^{y,y-N) 



d_ 
'd'z 



F{y,z;j, k) 



y-N 



F{y,y;j,k) — 



V.iy,z) + F{y,y-N;j,k) — 



V.{y,z) ) . 



-N 



It is easily verified that in this expression the first term 
vanishes identically, and after some calculation exploiting 
relations like 

{2iyN - zf 



2c2 



/(^, v) = {2vN - z) exp 

= f{z + 2N,y+l) , 

we find 

^^2^ VS^^^ = e-«'^' C{k) + e-''^^"^) C{N - k) 

(7.20) 



C(fc)=2iV^ [41.2/0(0,1^) e-^'^ 

-{2v-\f f^{N,v) e-t{^-k) 
-Ate? J2 ^ [fo{N, e-«"(^-'=) - 2/o(fc, i^) 

V 

-Atc^ dyVuh{y,v)e-^'\^-y\ , (7.21) 

•^0 
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where 



/o(y,!^) = exp 



{2vN - yf 



2c2 



(7.22) 



We now turn to the second part of the inhomogene- 
ity in Eq. (7.17), and we use Eq. (6.12) together with 
dc/dt = pilpo/{Nc) (cf. Eq. (7.2)), to write 



27rc-' 



N 



exp 



Thus 



27r e 



N 



N 



djoV*{jo,t\0) 



P^sPo Jo 

.^(,-f\3o-k\-f] ^ ^-f\N-jo-k\-q^{N-j)^ 



C{k) + e-«'(^-^) C{N - k) , 



(7.23) 



with 

C{k) = 4 ^ / djo 



N 



1 - 



• exp 



{2iyN+jof 
2c2 



(2i/iV+jo)' 



Q^\jo - k\ 



A short calculation shows that indeed 

C{k) = -C{k) . 



(7.24) 



(7.25) 



Thus the inhomogeneity in Eq. (7.17) vanishes and 
S{q,t;j,k,N) obeys the diffusion equation (7.18). 

Checking the boundary conditions (7.11) is an even 
simpler task. Direct calculation yields 



d_ 



S^^\q,t;j,k,N) 



r-N py 

= dy Vp^(y,2) 

Jo Jy-N ^ 

= fS^^Hq,t;0,k,N) 



d_ 
dj 



F{y,z;j, k) 



(7.26) 



d_ 
dj 



: [q, t; j, k,N) = -fS^''^ {q, t;N,k,N) , (7.27) 



and differentiating the integral equation (7.12), we find 

S'{q,t;j,k,N)=-^^S{q,t;j,k,N) 

d /■* /"^ 

= TP S^^'>{q,t;j,k,N) + dto djo V*{joM I 0) 
03 Jo Jo 

.|e-9^bo-fe| s'{q,t-to;j,0,N) (7.28) 
_e-f\^-jo-k\ S'{q,t-to,N-j,0,N)} . 



Now writing 

S'{q,t;0,k,N)=q^Si{q,t,k,N) , 
S'{q, t; N, k, N) = -q''S2{q, t, k, N) , (7.29) 

we find from Eqs. (7.26) - (7.28) 

Si{q,t,k,N) = S(-^\q,t;0,k,N) 

+ ldtoJ^ djo V*{jo,to I 0) {e-«"l^''-'=l 5i(g,i-io;0,7V) 
^^-f\N-oo-k\ S2{q,t-to,Q,N)] 

S2iq,t,k,N) = S^^\q,t;N,k,N) 

+ jdto djo V*{jo,to I 0) |e-9"IJ°-'=l S2{q,t-to;0,N) 
+g-f\N-jo-k\ Si^q^t-to,0,N)} . (7.30) 

This is exactly the system of equations obeyed by 
S{q,t;0,k,N) and S{q,t; N, k, N) (cf. Eq. (6.3)). The 
uniqueness of the solution together with Eq. (7.29) 
thus guarantees that the boundary conditions (7.11) are 
obeyed. 

We thus have shown that in the limit N 00 with 
t/Ts and q^R^ fixed, our theory reproduces the results of 
Doi and Edwards. 



VIII. NUMERICAL EVALUATION AND 
COMPARISON TO MONTE CARLO DATA 



A. Technicalities of solving the integral equations 

As shown by relations like Eqs. (6.19), (6.20), the 
natural measure of time in our theory is the parame- 
ter c = c{t). It measures the motion of the chain ends, 
i.e., the time dependence of tube destruction, and is de- 
fined by Eqs. (4.25), (3.16). In evaluating the theory, we 
therefore replace time by the variable 



X = X{t) = 



cit,N) 
N 



1 

N 



(8.1) 



Using Eq. (6.15), we write the integral equation (6.6) in 
the form 



Se{x) = S'^Pix) 



(8.2) 



^ dx' f^.^ , 1 



+2 — Ke [q'Nx', -j Se{X{t{x) - t{x'))) , 

where t = t{x) is the inverse function to a; = X{t). Se{x) 
and (x) denote the scattering fimctions 5*^; and , 
normalized with the static coherent structure function, 
e.g. 
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Se{x) 



SE{q,t;M,N) 
Sc{q,t = 0;M,M) 



(8.3) 



With corresponding notation, Eqs. (6.4) and (6.16) yield 
for the normahzed coherent structure function 



4 r dx' 



(p- Jo 



1 M 



(8.4) 

Se{X{t{x)-t{x'))) . 



In Eqs. (8.2), (8.4), we then transform from variables x' 
to x — X{t(x) — t{x')) to find equations of standard 
Volterra form, which are solved by discretizing x, x and 
iteration. We note that both S\< {x) and Sc {x) for 
X > 2 are negligibly small, less than 10~^, to be compared 
to the normalization 5c (0) = 1. Also the kernels JCe, 
exceed 10~^ only in the interval 0.1 < x' < 2.2. This 
allows for an accurate evaluation, simply using computer 
algebra. The numerical uncertainty of our final results 
is less than 0.5 %. In all our analysis, we used the same 
parameter values as in Ref. 0. Specifically, I^po = 1.23, 
p = 1/5, and = 2.364. (We recall that the precise 
values of p and ^s in fact are irrelevant for our numerical 
results.) 



B. Typical results 

We first want to illustrate the magnitude of the dif- 
ferent contributions to the normalized structure function 
Sc- Fig. 6 shows results for two very different values of 
wave number: q^R^g ~ 0.27 in Fig. 6 a, and q^R^g ~ 53 
in Fig. 6 b. The thick lines give the full results for Sc, 
including end effects and tube destruction. Long dashes 
represent Sc , i.e., the contribution without complete 
tube destruction. Short dashes represent the (normal- 
ized) contribution Si (Eq. (4.33)), which omits all end- 
effects and treats the chain as embedded in an infinitely 
long tube. The arrows point to the internal equilibration 
time fa, defined by Eq. (5.26): 



{N + lf 



Finally, the heavy slashes in the time axes give the repta- 
tion time defined as the first moment of the time depen- 
dent probability density of complete tube destruction: 



dto to 



(8.5) 



Here we use Eq. (6.19) for dV^'^yd{to). For long chains, 
the thus defined reptation time tends to the value given 
in Eq. (7.4). 

Fig. 6 a shows the typical behavior of Sc{q,t,N) for 
wave numbers which are too small to resolve the internal 

~(T) 

structure of the tube. Sc stays close to 1 up to times of 
order T3 and then rapidly drops to zero. The contribution 



of processes with complete tube destruction (dot-dashed 
line in Fig. 6 a) is very important, and end effects be- 
come visible only at a time where also tube destruction 
plays a role. Note that according to Eq. (6.4), the total 
structure function Sc is a sum of two independently cal- 
culated terms: Sc (long dashed) and the contribution 
of complete tube destruction (dot-dashed). These terms 
add up to a smooth curve (thick solid), an observation 
which demonstrates the consistency of our approach on 
the quantitative level. It is only for very short chains, 

^ 30, that these two contributions do not quite match. 
We trace this back to our approximate calculation of the 
kernel ICc- For such short chains, tube length fluctua- 
tions and internal relaxation presumably play a role also 
for the kernels. 

The limit of large wave vectors is illustrated with 
Fig. 6 b. Here configurations where the original tube 
has been destroyed, essentially do not contribute to the 

- ~(T) 

scattering. Indeed, in Fig. 6 b, the curves for Sc and Sc 
fall right on top of each other. However, end effects like 
tube length fluctuations have a strong influence, as shown 
by the deviations among the full line and the dashed line 
representing Si. They lead to a gradual decrease of S^ 
starting long before complete tube destruction becomes 
effective. 

In the previous section, we have shown that our the- 
ory asymptotically reduces to the primitive chain model 
of Doi and Edwards. To test the range of validity of the 
asymptotic result, we have evaluated our theory for the 
fairly large value N — 637 of the chain length (corre- 
sponding to a Monte Carlo chain of 640 beads, cf. Ref. 
0, Sect. II. C). Fig. 7 shows the results for the normal- 
ized coherent structure function Sc as function of logj^Q t 
for a set of wave vectors: q^R^ = 0.1, 1.0, and 10. The 
dashed lines give the asymptotic result (2.5), (2.6), where 
we used Eq. (7.8) for r^. Obviously, the time scales do 
not quite match: even for this long chain, the reptation 
time does not yet follow the iV'^-law. A shift of logj^Q t by 
—0.1, equivalent to a decrease of the time scale by 20 %, 
for small q, such that q^R^ ^ 1, brings the asymptotic re- 
sults close to the results of our full model. For q^Rg = 10, 
however, even with such a shift, there remains a definite 
difference: the result of the full theory initially decreases 
faster and approaches the shifted asymptotic curve only 
for i > 10^22- This is an effect of internal relaxation and 
tube length fluctuations. The absence of a visible mis- 
match in the shape of the curves for smaller 5- values just 
implies that with such small values again the structure 
of the tube cannot be resolved. 

To examine more closely the influence of internal re- 
laxation and tube length fluctuation, we in Figs. 8 and 
9 show results for the scattering from internal pieces of 
a chain. Fig. 8 shows results for g = 0.5 and a subchain 
of about 80 beads, corresponding to q^R^ ~ 3.29. The 
subchain is embedded as central piece in chains of dif- 
ferent lengths, precisely: (A^ = 77, total chain) , (A^ — 
157, M = 79) , {N = 317, M = 79) , (A^ = 00, M = 79) . 
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This figure illustrates the suppression of end-efi^ects with 
increasing N. The asymptotic result A'^ = oo is due to 
internal relaxation only. For finite N, the curves start 
to deviate from the asymptotic form as soon as wiggles 
created at a chain end have a nonnegligible probability to 
reach the central piece. The characteristic time for this 
process scales with {N — M)^. Fig. 9 shows results for 
q = 0.5, N = 317, and central pieces of lengths M = 317 
to M = 39. Due to tube length fluctuations, the nor- 
malized scattering function of the total chain initially 
decreases faster than the result for M = 159. Tube de- 
struction on the average reaches the sub-chains at times 
between i « 10^-^ {M = 159) or i « lO^-^ (M = 39), 
so that for large time regimes, the results for the sub- 
chains are not affected by tube length fluctuations. The 
decrease of the normalized scattering intensity with de- 
creasing length of the subchain rather is due to the fact 
that a shorter subchain leaves its initial position in the 
tube earlier, i.e., it is due to internal relaxation. 

In Figs. 8 and 9, we included results from a simula- 
tion of the Evans-Edwards model ||l^ . This model takes 
the chain configuration as a random walk on a cubic 
lattice and allows only moves of 'hairpin'-configurations 
rj_|_i — Tj = Vj-i — Tj as internal motion. An illustration 
for a two-dimensional system is shown in Fig. 1. We used 
the same implementation of the model as in our previ- 
ous work 1^, to which we refer for details. In comparing 
theory and simulations therefore all parameters are fixed 
by our previous analysis of segment motion. Since, how- 
ever, the new simulations lead to better statistics, we 
allow for some readjustment of the relation among t and 
the Monte Carlo time scale: t = 6.8 • 10~^ tMC instead 
of f = 6.09 • 10~^ tMC taken previously. This yields a 
shift of —0.048 of the logarithmic time scale and leaves 
our previous results essentially unaffected. 

As shown in Fig. 8, our theory in all details reproduces 
the time variations of the data, but the data systemat- 
ically lie somewhat below the theoretical curves. This 
is not due to our approximations which essentially only 
concern the treatment of tube length fluctuations. Con- 
sidering for instance the data for N — 317, we note that 
the deviations from the theory are strongest for t ^ 10^, 
where tube destruction and tube length fluctuations are 
irrelevant and our theory for the internal part M = 79 
is an exact evaluation of the reptation model. Further- 
more, the deviations are fairly independent of the lengths 
of the end pieces {N — M)/2. This suggests that we see 
some (non-universal) relaxation of the micro-structure. 
Clearly particle hopping, which is the elementary dynam- 
ics of the reptation model, is no faithful representation 
of the Monte Carlo hairpin dynamics on the microscopic 
level. The wave- vector |q| — 0.5 is large enough to resolve 
such details. Since the dynamic effects of micro-structure 
should saturate at larger times, this suggests that we 
should scale down the theoretical curves by some factor 
Bji < 1. This was done in Fig. 9, where Bj^ ranges from 
0.981 to 0.990, depending on M. For i > 10^, theory and 
data agree excellently. The same level of agreement can 



be reached for the data of Fig. 8. We therefore believe 
that our theory adequately describes the universal part 
of the coherent scattering function, including tube length 
fluctuations and (universal) internal relaxation. 

A more extensive presentation of simulation results, 
comparing with the present and previous theories, will 
be published elsewhere. Here we only note that we have 
taken data for values of |q| ranging from 1.0 to 0.1, and 
the g-dependence found for the initial deviation among 
theory and simulations strongly supports the interpreta- 
tion as micro-structure effects. 



IX. SUMMARY AND CONCLUSIONS 

In this work, we have exploited the pure repta- 
tion model to calculate the coherent structure function 
Sc{q,t; M, N) of a flexible chain moving through an ar- 
ray of impenetrable topological obstacles. Our analysis 
is rigorous for a subchain (length M) in the interior of 
an infinitely long chain {N oo). This allows for a 
detailed comparison with an approach where the inte- 
rior motion of the chain is modeled as one dimensional 
Rouse motion along a coiled tube. We found that the lat- 
ter model starts from unphysical non-equilibrium initial 
conditions, which relax only on the scale of the Rouse 
time T2 (M) of the subchain. This relaxation completely 
distorts the time dependence of Sc{q, t; M, N — 00) for 
times t<T2{M). Only for times t > T2(M), this model 
is equivalent to the reptation model. If applied to the 
total chain, 'local relaxation' calculated as Rouse motion 
in a tube therefore is unphysical. A realistic system may 
show some relaxation which is specific to the microscopic 
dynamics, and which is not contained in the pure repta- 
tion model. However, our analysis sheds strong doubts 
on an interpretation of such non-universal effects within 
the framework of the model of a Rouse chain in a tube. 

To evaluate the total structure function for all times, 
we have derived integral equations which split Sc into 

(T) 

a contribution Sc of configurations where some part of 
the initial tube still exists, and the remainder. The kernel 
and in particular the inhomogeneities (like ) of these 
equations cannot be calculated rigorously. They involve 
distribution functions coupling the motion of a given seg- 
ment to tube renewal, which is a non-Markovian process 
with memory time of the order of the Rouse time T2{N). 
To calculate the functional form of these distributions, 
we used a random walk approximation. We thus at each 
instant of time replaced the correlated process by an un- 
corrclated process which as closely as possible reproduces 
the instantaneous distributions of the correlated process. 
This 'mean hopping rate' approximation introduces func- 
tions c = c(t), a = a{j, t), and b = b{j, t) which appear as 
parameters in the distribution functions and can be cal- 
culated from the microscopic hopping process of spared 
length in the reptation model. They also have a sim- 
ple physical meaning: c{i) measures the average extent 
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of tube destruction, a{j, t) describes the couphng of mo- 
tion of segment j to tube destruction, and t) takes 
care of the inhomogeneity of the effective segment mo- 
bihty along the chain which arises from the fact that the 
mobile units, i.e., the wiggles of spared length, can be 
created and destroyed only at the chain ends. 

In the limit of long chains and time t ^ T2, the param- 
eters a and b tend to 1 irrespective of j, and c ~ t^^^. Our 
theory then reproduces the results of the 'primitive chain' 
model. Our proof of this asymptotic result amounts to 
a derivation of the primitive chain model from micro- 
scopic reptation dynamics. Combined with the analysis 
of the internal motion for t <C T2, this result allows for a 
mapping of the microscopic parameters of our reptation 
model to the more commonly used Rouse-type parame- 
ters. 

The parameter functions a{j,t), b{j,t), and c{t) ap- 
proach their asymptotic behavior only slowly, and it 
needs chain lengths of order N/Ng > 300 to find a time 
region where the primitive chain model is valid. In par- 
ticular, for shorter chains the reptation time, if extracted 
by fitting the Doi-Edwards result for the primitive chain 
model to the large time behavior of Sc{q, t, N) in the rep- 
tation model, does not obey the asymptotic power law 
T3 N^. As will be shown in Ref. Q, it rather exhibits 
the well known behavior T3 ~ N^"^" , with an effective 
exponent ZeS > 3. The deviation of a(j, t), b{j,t), and 



c{t) from their asymptotic primitive chain behavior incor- 
porates the effect of internal relaxation and tube length 
fluctuations. Our numerical evaluation of the full theory 
illustrates that these effects in general are quite impor- 
tant. In particular, we find a clear difference among the 
time variation of scattering from the total chain com- 
pared to scattering from internal pieces. The latter are 
less influenced by tube length fluctuations but are more 
strongly affected by internal relaxation. This leads to a 
peculiar behavior of Sc{q,t;M,N) with varying length 
M of the internal piece, as shown in Fig. 9. Quite gen- 
erally, for the total chain {M = N) it is the tube length 
fluctuations, that determine Sc for times up to < ~ 10 

All our quantitative numerical results are well sup- 
ported by simulations of pure reptation, exploiting the 
lattice model of Evans and Edwards. In view of the 
unavoidable approximations inherent in the theory, the 
quantitative agreement is quite remarkable. It suggests 
that our mean hopping rate approximation adequately 
takes care of the coupling among internal relaxation, 
tube length fluctuations, and global creep. A more ex- 
tensive comparison to Monte Carlo data including a nu- 
merical parametrization of our analytical results will be 
published elsewhere 0. 
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APPENDIX A: RANDOM WALK MODEL FOR DISTRIBUTION FUNCTIONS 



1. The function ^^'^^(nm, no; 



As explained in Sect. IV. C, we consider a random walk n'{s) on the integer numbers, with hopping rate p'. The 
walk starts at n'(0) = and ends at n'{t) = no- It is restricted to the interval [rim — N' + l,nm], with absorbing 
boundary conditions. To simplify the notation, we shift the interval by N' — rim to 



J = [1,7V'] 



(Al) 



with the starting point of the walk h'q = N' — Um, and the endpoint n'{t) = no + N' — nm- The hopping matrix of 
the walk takes the form 



WiN'hf = (1 - 2p')%' +p'{Sj,j'+i + 5j^,j-i), e I . 



(A2) 



'^n^ o('^"i) no; t) gives the weight of the walk under the constraint that n'(s) = nm is attained for at least one s G [0, t]. 
It is easily found as 



|n| + n 



e 



n — \n\ 



nm + N' -I 



{(W^*(iV'))l+„^-„,l+„^ - (1 - Sn^fi) {W\N' - l))„^-„,„^} . 



W^{N') can be written as 



N' 



N' 



UK] , TTKJ 

sm — sm ■ 



N' + l N' + 1 



1 - 4p' sin^ 



2{N' + 1) 



(A3) 



(A4) 



For 1, the last factor can be replaced by exp (^—p'tj^T-fjyi^ , and a little calculation yields 



|n| + n 



e 



n — \n\ 



-nm + N' -1 



■' + 1 ^ 



N' + 



TTK n TTK (2nm + 2 — n) 

COS — — - — cos - 



N' + 1 



N' + l 



1 ^ 



1^ 7T__K_ 
N'2 



K = l 



TTK n TTK (2nm — n) 

COS —-r. cos - 



N' 



N' 



(A5) 



In extending the sums to infinity, we neglect terms of order exp(— Tr^p't). 

(T) 

Now we note that 'PmaxO rapidly decreases for t > T3. As mentioned in Sect. III.B, the effective hopping rate for 
i ~ T3 » T2 behaves as p' ~ 1/N ~ l/N', and as a consequence, the argument of the exponential in Eq. (A. 5) takes 
the form —const ti^t/T^. Thus V^^o yields a relevant contribution only for smaller times: t < T3, and for treating 
this time range, it is preferable to replace the summations in Eq. (A. 5) by their Poisson transform. This yields 



'P^Loi^m, n;t) = e[nm + 



\n\ + n 



e 



n — \n\ 



nm + N' -I 



y/Airp't 



■5 »p(-i^Knr)h(-|K-f)-^)-' 



(A6) 



exp 



1 / .r/ ( 2j^ + 2/ n\ [y + \f\ II 

(^^ + - 2 ) j h ("^ ('^^ + - 2 ) - j - I • 



Now the Gaussian prefactors of the square brackets allow for an essential contribution only for p't > iV'^, and in this 
region, the square brackets can be replaced by the linear approximation 
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exp 
2v 



2v 



vN' 



2u + 2 



vN' + exp ; — vN' + ; 

p't\ 2) ^ [ p'i V 

2(z^ + l) / AW n\ 



2j 



p't 



p't 



Eq. (4.20) from Sect. IV. C is the result, which is correct to leading order in \/N' . The neglect of l/A'^'-corrections is 
consistent with treating segment indices as continuous. 



2. The function V*{jo,to \ 0) 

By definition, V*{jo, to \ 0) gives the probability that the initial tube is destroyed completely at time step to, with 
jo being the last point, occupied by chain end (see Sect. VI. A). In the random walk model, this probability is given 
by the weight of a walk starting at n'(0) = and ending at n'{to) = jo = jo/^s- The point is reached at t = to 
for the first time, but the point jg — N' + 1 is attained for some s G [0,to — 1]. The walk is restricted to the interval 
[jo — N' + l,j'o]. Shifting the interval by N' — j^, we can express this weight as 



V*ijo,to \0)=p' UW^^-^N' - l))jv'-i,jv'-j' - {W'°-\N' - 2))jv'-2 



,N'-i'-l 



(A7) 



where the prefactor p' gives the probability of the last step, leading from A'"' — 1 to N' (in the shifted walk). Using 
the explicit expression (A. 4) for W* and exploiting i = to — 1 1, we find 
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N' - 1 
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(A8) 



correct up to exponentially small terms (cf. Eq. (A. 5)). The Poisson transform yields 
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(A9) 



As in Eq. (A. 6), we can expand in the square brackets, keeping the first non- vanishing terms which here are of second 
order. Identifying now N' = N/Ig and j^ = jo/^s, Eq. (6.11) in Sect. VI. C results. 



APPENDIX B: EXPRESSION FOR vf \nj;t) 

According to Eqs. (4.28) and (4.30), Vj (% ; t) can be determined by integrating Pj^g^j over Um (or y, equivalently) . 
This integral can be carried through analytically to yield 
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- (2 - a^)v - - erfc 



(2-a2)zy+- erfc 



az + a2vN + 



az + a2vN - N 



+ ua{z - 2uaN) 



2{az + a2yN) erfc 
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[az + a2yN + N) erfc 



- {az + a2vN - N) erfc 



az + a2vN + N 



/02 

az + a2vN - iv' 



— 1^0(2; — 2vaN) 



/ (az + a2J^Af)M / {az + a2i'N + Nf 
2 exp — exp 



02 



02 



- exp 



{az + a2uN - Nf 
a2 



(Bl) 



with 02 = 2(1 - a^). 



APPENDIX C: MODELING A ROUSE CHAIN IN A COILED TUBE 

The model describes the internal dynamics of the reptating chain as that of a one-dimensional Rouse chain, stretched 
so as to span the contour length of the tube. The potential energy takes the form 



772Zl(2;j-^J-i)^-7(^JV-a;o) , (CI) 



where the Xj, j — 0, . . . ,N are the bead positions, i measures the mean segment size, and h/i is the stretching force 
acting on the end beads. The average extension of the chain is easily calculated as 

L = {xn -xq) = 2ehN . (C2) 

The dynamics of the chain is given by a Langevin equation. 

d_ d V 

It ~ dx] 

The fluctuating force is Gaussian distributed: 



^.=-70 — 17^+?, (C3) 



P[^j{t)]=M-' exp 



1 /■+°° 



(C4) 



It is a standard exercise to calculate the dynamical structure functions. Indeed, for this system of coupled harmonic 
oscillators, the stretching force does not influence the dynamics but changes the static prefactor only. For the scattering 
from a pair {j, k) of beads, one finds 



where 

Djk{t) = ^ + e'\j-k\ 
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Using the approximate form of w^, we neglect some exponentially small micro-structure effects. 

In the analysis of Sect. V, we need this result for two segments deep inside a very long chain, for times small 
compared to the Rouse relaxation time of the total chain. Writing j = N/2 + j, k = N/2 + k and taking the limit 
N 00, with 7oV^^ ^ 1 fixed, one finds 



Djk{t)=f \j-k\+£^/2^tg 



|j-fcK 



(C8) 



where 



z erfc2; . 



(C9) 



Now assume that the chain is embedded in a coiled tube, consisting of N/N^ segments of fixed length It- The 
contour length of the tube equals the length of the stretched chain: 



N 

lT — =L = 2thN. 

N,, 



(CIO) 



The end-to-end distance of the tube i?g = must match the end-to-end distance of the physical chain. The 

potential energy of the latter is given by the three-dimensional version of Eq. (C.l) in the absence of the stretching 
force, which results in B? = 6£'^N. Thus 



(Cll) 



and Eq. (C.IO) yields 



2/1= W— . 

Ne. 



(C12) 



To calculate the scattering from segments j and k of the stretched one-dimensional Rouse chain embedded in the 
tube, we write 



Tube 



J — OO \ 

= r^dx /e^q(r.(*)-r.(0)) 

^ — OO 



Xk{Q) - x) 



Tube 



Xj{t)-Xk{0)=x 



5{xj{t) - Xk{Q) - x) 



(C13) 



Here the first factor is to be calculated under the constraint that the points Yj{t) and rfc(O) have distance x measured 
along the tube. It is thus given by the static correlation function of a chain oIx/It segments of fixed length It: 



,iq(r,(t)-rfc(0))^ 



(C14) 



The second factor in Eq. (C.13) is the (1-dimcnsional) Fourier transform of S'^^'^^ (Eq. (C.5)) 



(5(a;j(t)-a;fe(0)-a;) ) = 



+ 00 



dp 



= (47rD,fc(t))-i/2exp 



{x - 2he{j - k)f 



Substituting Eqs. (C.14) and (C.15) into Eq. (C.13), we can carry out the integral to find 



(C15) 
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S^^\q,t;j,k) = i e<5'{e2^^ erfc (g + A) + e'^AQ erfc(Q- A)} 



(C16) 



-Djkit) 



(C17) 



A = 



(C18) 



If we take for Djk{t) the result (C.8), the variables Q and A reduce to Q and A given in Eq. (5.15), and the result 
(C.16) becomes identical to the expression (5.14). 

A final remark on De Gennes' approximation may be appropriate. The derivation starts from Eq. (C.5) with 
Djk taken from Eq. (C.8). Aiming directly at the coherent structure function, one integrates this expression over j 
and k. If we ignore end effects, this yields 



1 / djdk S^^^\p,t-j,k) 



ds exp 



where s = \j — k\. To evaluate this integral analytically, one expands exp —p'^i y/2jot g (^--^|=j up to first order. 

The remaining steps closely follow our derivation given above and result in the form (2.12) of the 'local' contribution 
to the coherent structure function. It should be noted that the expansion is valid only for \/2rfot <C 1. The analysis 
supposes that t is small compared to the Rouse relaxation time, and in this time regime, ^/2joi is of the order of the 

mean square distance moved by a segment along the tube. Thus the condition p'^l ^/2jot k p^ {x.j [t) — Xj (0))^ ^ 'ti I 

implies that the wave number p cannot resolve the motion of a segment. 
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FIG. 1. A realization of reptational dynamics in the Evans- Edwards lattice model (two dimensional illustration). The crosses 
denote impenetrable obstacles which allow only for 'hairpin' moves as shown by the sequence of pictures. The hairpins represent 
the wiggles of spared length. 
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b) 

FIG. 2. Schematic drawings of the cases discussed in the text. The full line represents the unfolded initial tube. Broken 
lines represent the new end pieces of the unfolded tube created up to time t. 



30 



logioit) 



FIG. 3. 




, 

Coherent normalized structure function Sc{q,t; M, oo) of the central piece of an infinitely long chain for 
q~Ug(M) — 50, M — 200. Fat line: reptation result (Eq. (5.4)); long dashes and short dashes: Rouse chain in a coiled 
tube. For the latter curve, a constant has been subtracted such that Sc{q, 0; M, oo) = 1. Medium size dashes: De Gennes' 
approximation (Eq. (5.11)). The arrow points to T2{M = 200). 




FIG. 4. Initial (0) and final (t) configuration of the chain (fat lines), together with the configuration at time to (thin line). 
At time to, the chain leaves the last piece of the initial tube, with one chain end at the position of bead jo in the initial tube. 
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FIG. 5. Distribution functions for complete tube destruction, a) Probability density of complete tube destruction as function 

of 2/Af2 - t/Tz. Normalization: d{2/N^) P*(0, iV, iV) = 1. 

b) Probability density of complete tube destruction as function of the position jo/N of the final segment of the original tube: 
d{2/N^) V{jo/N) = dto 'P*{jo,to \ 0). The chain leaves the tube with end 0. The values of chosen are indicated by 

arrows in Fig. 5 a. 
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FIG. 6. Normalized coherent scattering function in different wave number regions as function of logj„(f). a) — 0.01, 
N = M — 157; b) = 1.0, N = M = 317. Tfie thick solid lines give the full functions Be- Long dashes represent S'^\ which 
in b) coincides with Sc- Short dashes are the results neglecting all end effects. The dot-dashed line in a) is the contribution of 
complete tube destruction. See the text for further explanations. 
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1 2 3 4 5 6 log^Q(t) 

FIG. 9. Results for Sc, q = 0.5, N = 317. Theoretical results (solid lines) for central sub-chains of lengths M =39, 79, and 
159 are compared to simulations (dots). Results for the total chain (M = 317) are also shown. 
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